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Abstract 



In this more pedagogical study we want to elucidate on stochastic aspects 
inherent to the (non-)equilibrium real time Green's function description (or 
'closed time path Green's function' - CTPGF) of transport equations, the 
so called 'Kadanoff-Baym equations'. As a toy model we couple a free scalar 
boson quantum field to an exemplaric heat bath with some given temperature 
T. It will be shown in detail that the emerging transport equations have to 
be understood as the ensemble average over stochastic equations of Langevin 
type. This corresponds to the equivalence of the influence functional approach 
by Feynman and Vernon and the CTP technique. The former, however, gives 
a more intuitive physical picture. In particular the physical role of (quantum) 
noise and the connection of its correlation kernel to the Kadanoff-Baym equa- 
tions will be discussed. The inherent presence of noise and dissipation related 
by the fluctuation-dissipation theorem guarantees that the modes or particles 
become thermally populated on average in the long-time limit. For long wave- 
length modes with momenta \k\ <C T the emerging wave equation do behave 
nearly as classical. On the other hand, a kinetic transport description can 
be obtained in the semi-classical particle regime. Including fluctuations, its 
form resembles that of a phenomenological Boltzmann-Langevin description. 
However, we will point out some severe discrepancies in comparison to the 
Boltzmann- Langevin scheme. As a further byproduct we also note how the 
occurrence of so called pinch singularities is circumvented by a clear physical 
necessity of damping within the one-particle propagator. 

PACS numbers: 11.10. Wx, ll.15.Kc, 5.40.+j, 5.60.+W 
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I. INTRODUCTION AND MOTIVATION 



Non-equilibrium many body theory or quantum field theory has become a major topic of 
research for describing various transport processes in nuclear physics, in cosmological particle 
physics or more generally in quantum dissipative systems. A very powerful diagrammatic 
tool is given by the 'Schwinger-Keldysh' |l|-|3| or 'closed time path' (CTP) technique by 
means of non-equilibrium Green's functions for describing a quantum system also beyond 
thermal equilibrium M. For an equilibrium situation this technique is equivalent to the 
real time description of finite temperature field theory PJ^,§|- The resulting equations of 
motion (by various approximations), the so called Kadanoff-Baym equations have to 
be considered as an ensemble average over the initial density matrix p^(t ) = \i)pfj(j\ 
characterizing the preparation of the initial state of the system. Typically, if the system 
is close to thermal equilibrium, the (initial or resulting) density matrix allows for thermal 
fluctuations which should be inherent to the transport process under consideration. If fluc- 
tuations are physically present in the course of the evolution, then, according to the famous 
fluctuation-dissipation theorem, also dissipation must be present. Or, in turn, if the (quan- 
tum) system behaves dissipatively, as a consequence, there must exist fluctuations. Again, 
the Kadanoff-Baym equations have to be understood as an ensemble averaging over all the 
possible fluctuations. Accordingly, there should be a stochastic component inherent to these 
transport processes. 

This inherent stochastic aspect of the Kadanoff-Baym equations is what we want to point 
out in detail below. An important ingredient to these equations and thus to the one-particle 
propagator will consist of the ensemble average over 'noise' fluctuations. 

It might well be that some or most of our conclusions are already known and can be found 
in various studies. However, we feel that there exists no clear and thorough study which 
contains all the points we will present in the following. We believe that our investigation 
can give some new physical insight on the structure of the Kadanoff-Baym equations. 

As a primary intention and motivation for our further study let us make a seemingly 
jump by reviewing the description of (classical) Brownian motion. (The intimate connection 
of 'quantum Brownian motion' to the structure of the Kadanoff-Baym equations will be the 
basic outcome.) A heavy particle with mass M is placed in a thermal environment of fast 
moving particles (being much lighter, i.e. their masses m -C M). The equation of motion 
of the heavy particle can then phenomenologically be described by means of a standard 
Langevin equation: 



Here 7 describes the (Stokes) damping coefficient resulting from the average interaction 
with the lighter particles. £(t) is the random (Langevin) force (or 'noise') exerted by the 
surrounding medium (the 'heat bath'), the external agencies. Since this random force £(t) 
represents instantaneous effects of interactions (i.e. collisions) with many particles, it is 
plausible that is a stochastic force with a Gaussian distribution having an average of 
zero. 

It seems worthwhile to point out the specific nature of such a stochastic equation of mo- 
tion: The stochastic force term makes the trajectory x(t) a stochastic quantity as well. 



Mx + 7X = £(t) . 
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Thus, to any possible ensemble {£\(i),t = — > 00} we get a corresponding ensemble 
{x\(t),t = — > 00}. The physical description of the system x(t) results from piling up 
such trajectories. This means that we have to deal with ensembles of stochastic trajec- 
tories {x\(t)} instead of one trajectory x(t), as it would be the case for the solution of a 
deterministic equation of motion. 

Further, one typically assumes that possesses an infinitely short time correlation 
(which physically corresponds to a clear separation among the time scales of relaxation of 
the heavy particle and the light particles), i.e. that it is a Markovian process and the noise 
being white. £(t) is thus completely characterized by the following moments 

«*(*)» = 

(mm)) = m-t') . (1.2) 

((...)) labels the average over all possible realizations of the stochastic variable £(£). All 
higher moments follow as a factorization of 'truncated' second moments and a sum over all 
possible permutations from the property of the (assumed) Gaussian distribution. 

The friction coefficient 7 and the random force have the same microscopic origin. 
In fact, it turns out that there is a relation between them, called the fluctuation- dissipation 
theorem (FDT) flOf . (Such a relation will explicitly shown to exist within the CTPGF 
approach.) As a motivation for its microscopic existence, one should recall that for long 
times the average energy of the heavy particle should become ((E kin )} = ((p 2 )} / (2M) = T/2. 

The explicit solution of ( |1.1| ) is readily found to read (/3 := 7/M) 

t 

pit) = J dt' e-^-^Ot') + e-^p(O) , (1.3) 


so that with ( |1.2|) one has 

((P 2 (t))) = ^ (l _ e-^) + e-^«p 2 (0)» L (1.4) 

and thus the desired connection follows as 

I = 2M/3T = 2 7 T . (1.5) 

The strength of the noise distribution / is related to the friction coefficient 7, which is a 
rather simple version of the FDT. (Such a relation is often been referred in the literature as 
the second FDT.) The microscopic derivation of the friction 7 and of the residual random 
force £(t) is a longstanding problem of statistical mechanics. For a rather clear presentation 
and derivation for a classical as well as quantum mechanical system we refer to the article 



by Cortes et al. fllT |. One can also convince oneself that the higher moments ((p n )) turn out 



to be the one for a free particle with a canonical distribution, i.e. 

«p n (*)» ^ (p n U = d 3 pe~^p n (1.6) 

with 

Zq := J d 3 pe~ 



v 1 

2MT 
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This follows immediately by making use of the Gaussian distribution and the resulting 
factorization of the higher moments. 

If there exists a clear separation among the time scales of the slow heavy particle and 
the fast light particles, a Markovian Langevin equation with white noise can be considered 
as being appropriate. However, this might not truly be the case for various applications like 
in nuclear physics [|T2"| , |T3ll . In our later field theoretical discussion the Markov assumption is 
too restrictive. We are thus forced to generalize the above Langevin equation ( |1.1[ ) to allow 
for finite memory, i.e. retarded friction. It reads 

t 

Mx + 2 J dt'V(t-t')x(t') = f(t) . (1.7) 

— oo 

(The Markov approximation is thus obtained if F(t — t') = 7<5(t — £').) We have shifted 
the lower bound of the integration to — oo which is valid for the discussion of the long-time 
behavior. Tacitly we assume that the friction kernel T(t) = T(—t) is 'well behaved' (meaning 
that its Fourier transform T(u) > 0) and has only some finite extension in time. On the 
other hand, in the (semi-)classical regime fill the time correlation of the stochastic force 



follows as 

i( t -t') : = (mm)) = 2rr(f-o (= amm))) ■ o~*) 

A noise correlation being non-local (i.e. not white) is called colored noise. 

To see that this relation is still in accordance with the equipartition condition 
((p 2 ))/(2M) = T/2, we define the retarded Green's function 

2 } 

G TCt (t) = -— J dt'T(t-t')G vct (t') + 6{t) 

— oo 

with the boundary condition: G ret (t < 0) = 

-> GretM = 7 WTT\ T (1-9) 

and the advanced Green's function (G rct (t) = G m {— t)) 



oo 



-Gav(t) = -— J dt'T(t-t')G av (t') + 5(t) 
t 

with the boundary condition: G av (t > 0) = 



- = 7-— — ;i ufA — v = ■ (i.io) 

(Note: Due to the symmetry condition in ( |1.8| ) it is f (a>) = r(— u) = and we assume 

F(u;) > for 'true' dissipation. As indicated in (|1.9Q , (|1.10|) , an infinitesimal small width 
e > has to be added if f (oo) = for some frequencies ui so that the retarded propagator 
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only has poles in the lower complex half plane and the advanced propagator, respectively, 
only has poles in the upper complex half plane.) In the equations above P denotes the 
principal value of the integral and T is the Fourier transform of T. By means of ( |1.9| ), ( |1.10| ) 
one finds in the long-time limit the desired property 



du - 



TM j — 2 IMB. , = TM . (1.11) 

it M J i, j — l,y 



rM j u-u' J ' M 2 

(We have not given the proof for the last equality. This follows rigorously from Cauchy's 
theorem in the extended complex plane by making use of the analytic properties of the 
respective propagators and assuming T(u) < \uj\~ v , r\ > 0, for \u\ — > oo.) Hence, irrespective 
of the detailed form of the friction kernel T(t), the equipartition condition is automatically 
fulfilled as long as the noise kernel I(t) fulfills the relation ( |1.8| ). 

After this somewhat lengthy discussion of classical Brownian motion we now want to 
return to quantum field theory and already point out some similarities. One of the major 
topics in quantum field theory at finite temperature or near thermal equilibrium concerns 
the evolution and behavior of the long wavelength modes. These often lie entirely in the non- 
perturbative regime. Therefore solutions of the classical field equations in Minkowski space 
have been widely used in recent years to describe long-distance properties of quantum fields 
that require a non-perturbative analysis (for some applications see references cited in [14|]). 



A justification of the classical treatment of the long-distance dynamics of bosonic quantum 
fields at high temperature is based on the observation that the average thermal amplitude of 
low-momentum modes is large. For a weakly coupled quantum field the occupation number 
of a mode with wave vector p and frequency up is given by the Bose distribution 

n(up) = (f*/ T -l)~\ (1.12) 

For temperatures T much higher than the (dynamical) mass scale m* of the quantum field, 
the occupation number for the low-momentum modes \p\ -C T becomes large and approaches 
the classical equipartition limit 

n{ Uf ) » 1. (1.13) 

The classical field equations (which can be understood as a coherent state approximation 
to the full quantum field theory) should provide a good approximation for the dynamics 
of such highly occupied modes. However, the thermodynamics of a classical field can only 
be defined if a 3-dimensional ultraviolet cut-off k c is imposed on the momentum p such as 
a finite lattice spacing a. In addition, the long wavelength modes do still interact with 
the short wavelength modes (of the scale of the temperature). In a recent paper of one of 
us it was shown how to construct an effective semi-classical action for describing not 
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only the classical behavior of the long wavelength modes below some given cutoff k c , but 
taking into account also perturbatively the interaction among the soft and hard modes. By 
integrating out the 'influence' of the hard modes on the two-loop level for standard 4 - 
theory the emerging semi-classical equations of motion for the soft fields become stochastic 
equations of motion of Langevin type: 



□ 2 + m 2 th + S ret ® )4> soft = f(t) , (1.14) 



resembling the analogous expression to eqs. (|1.1|) and (|1.7f) . Splitting S ret into real and 



imaginary parts it follows that the (onshell) damping rate in the weak coupling regime 
(Ep 2> |E ret (p, E p )\) is identified as 

lm (Z Tet (p, E p )) 

lif) = 1 , (1-15) 

being responsible for dissipation. In addition, also the noise correlation function as the ana- 
logue to (1.8) follows microscopically We will see how this works in detail below. The 



hard modes thus act as a heat bath. They also guarantee that the soft modes become, on 
average, thermally populated with the same temperature as that of the heat bath. Equiv- 
alently to ( pL.ll|) we will see that a similar relation exists for the average fluctuation of the 
amplitude squared. For the semi-classical regime where \p\ <C T it will take the form 

((|0(p,t -> oo)| 2 » = V J ^D iet (p,u)I(p,uj)D av (p,uj) 

^^T^^ r N(p,t)^^ r (l + N(p,t)) , (1.16) 
E p Ep Ep 

where I(p, u) describes the Fourier transform of the time correlation function of the noise 
distribution and V is the volume of the system. Such kind of Langevin description for 
the propagation of long wavelength modes in non Abelian gauge theories has recently been 



proposed by Huet and Son |L5| to be used in the context of cosmological baryogenesis. 

The derivation of the effective stochastic equation of motion has been carried out by 
means of the influence functional approach of Feynman and Vernon [[H]]. The basic idea is 
to integrate out (perturbatively) to some order (typically quadratic) the heat bath degrees 
of freedom and thus obtain an effective action for the system degrees of freedom. One 
can explicitly see how dissipation and noise occurs and how the FDT emerges from a purely 
microscopic point of view. In fact, the idea of employing the Feynman and Vernon functional 
approach was used quite early to derive an effective Langevin equation for quantum Brownian 



motion |T8|-pO|. One can also formulate the effect of the heat bath on the evolution of the 
system degrees of freedom by means of the CTPGF technique. The main success of the 
CTPGF formalism lies in its 'usefulness' obeying strictly diagrammatic rules so that higher 
order contributions (i.e. diagrams) can, in principle, be straightforwardly calculated (in 
analogy to the Matsubara rules). In the CTPGF formalism, however, only the ensemble 
averaged characteristic of expectation values is obtained (as we will see). In fact, one can 
think of the individual stochastic trajectories being sampled and averaged over the influence 
of noise. This is what we want to discuss in depth. The equivalence of both approaches has 
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indeed been pointed out already in |21] on a quite formal level. However, we believe that 
our discussion is inspired on more physical intuition and thus hopefully provides some new 
insight on the dissipative and stochastic character of the CTPGF approach. 

In the following we discuss both approaches for a free scalar field theory interacting 
(perturbatively) with a heat bath. In the next section we review the CTPGF technique and 
evaluate the average characteristic properties of the CTP propagator resulting from the (non- 
equilibrium) equations of motion, the Kadanoff-Baym equations. We obtain the Feynman- 
Vernon approach by a rather simple rearrangement of the interaction kernels stemming 
from the heat bath. With this, of course, the equivalence has already been shown. We 
identify the dynamically generated mass shift, dissipation and fluctuation terms and show 
that the fluctuation-dissipation theorem emerges naturally and is thus stated in microscopic 
terms. In the third section the (quantum) noise distribution and its (in this case Gaussian) 
distribution will be introduced. With this at hand we analyze the 'noisy' propagator for a 
given noise sequence and then explicitly show that when averaging over the noise fluctuation, 
the CTP propagator results naturally. We thus can point out the inherent meaning of the 
noise correlation function on the propagator and on the CTP description. In section [TV] 
an attempt will be undertaken to derive an effective Boltzmann-Langevin equation for the 
phase-space occupation number of the bosonic particles. In this respect, we first show that 
the noise drives also the thermal fluctuations in the occupation number. We comment on the 
interpretational difficulties of the effective equation and on the quantum aspects of noise. As 
a byproduct of our study we also see how so called pinch singularities do emerge in a strictly 
perturbative sense and how this problem is cured naturally in this case (our approach for the 
system degrees of freedom is non-perturbative in the sense that the perturbative contribution 
(the 'influence') of the bath is resummed by a Dyson-Schwinger equation). This is discussed 
in section |V|. We close our findings with a summary. Four appendices are added for some 
technical issues. 



II. CLOSED TIME PATH GREEN'S FUNCTION TECHNIQUE 



We start with the closed time path action for a scalar field coupled to an exemplaric 
heat bath with temperature T, 



S 



1 r c 
, (fx 

2 

- d 4 x 



(f> c (x) (-□„ -m 2 )4> c (x) 



1 f c 



dfai d 4 x 2 [4>c{xi) S c (xi, x 2 ) 4> c (x 2 )] 



cj) + {x) {-U x - m 2 ) 4> + {x) - (f>~(x) (-n x - m 2 ) (j)"(x) 



-jd 4 x 1 d 4 x 2 + (xi) S ++ (xi,x 2 ) + (x 2 ) + + (xi) S + (x 1 ,x 2 )(fi (x 2 ) 

+ <j>~{xi) E~ + (X!,X 2 ) (j) + (x 2 ) + (f)~{xi) E~~(xi,X 2 ) 4>~(x 2 ) 



+ (-□ - m 2 )0+ - 0- 

- </>+ £++ <p + - </>+ 57 



-□ - m 2 ) <\T 

cj)- - </)- -IT+ <p + - (f>- E— <p~ 



(2.1) 



The interaction among the system and the heat bath is stated by an interaction kernel 
involving a self energy (or 'mass') operator E resulting effectively from integrating out the 
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heat bath degrees of freedom. Schematically this is sketched in fig. [j]. We assume that £ 
is solely determined by the properties of the heat bath. Ideally this should hold for the 
linear response regime. In appendix [A] we give a short derivation of eq. ( |2.1| ) resulting from 
a simple microscopic model within linear response. (We recommend the reader to start 
with appendix [AJ if (s)he is not too familiar with the real time Green's function technique.) 
In the first line of (|2.1|) we have stated the closed time path action within the Keldysh 
representation with the time integration along the closed time contour C (see also appendix 
[A]). In the second line we have disentangled the time integration on the contour to ordinary 
time integrals on the usual real time axis. This corresponds to the usual conversion from 
the closed time description to the 2 x 2-matrix notation |3|-[5],|7|] . The last equation defines 
obvious abbreviations we will use below. If not explicitly stated otherwise the integrations 
are defined as 

dSc := Idt ld 3 x (2.2) 



with to being some initial time at which the system starts to evolve from a specified ini- 
tial density matrix p( l \to). In (|2.1| ) the self energy contribution from the heat bath is 
parametrized by the four self energy parts 

Z ++ (xi, x 2 ) = e(tx - t 2 ) S > (x 1 , x 2 ) + Q(t 2 - t x ) S < (x 1 , x 2 ) = E n (xi,x 2 ) , (2.3) 

£~ (x ls x 2 ) = &(t 2 - h) S > (xi, x 2 ) + 0(ti - t 2 ) £<(£!, x 2 ) = E 22 (x 1? x 2 ) , (2.4) 

E + ~(xx,x 2 ) = —T* < (xx,x 2 ) = -S 21 (xi,x 2 ), (2.5) 

£~ + (xi,x 2 ) = -S > (x i ,x 2 ) = -S 12 (xi,x 2 ), (2.6) 

where tx and t 2 are the time components of the vectors xx and x 2 , respectively. On the 
right hand side of each of the equations we have stated the corresponding expressions of 
the four self energy parts when using the contour notation as employed in appendix [A] (see 
eqs. ( |A19| ). Without specifying S > and E < we can already observe here that 

E ++ + £ + £+- + 5T+ = (2.7) 

holds. 

The self energies S > and E < generally fulfill the following relation in case of a scalar 
boson field theory 

S > (xi,x 2 ) = £<(x 2 ,xi). (2.8) 

As an example one can easily convince oneself that this relation is fulfilled for the explicit 
expressions stated in eqs. ( |A19| ). Relation ( |2.8| ) is a special version (being true for a real 
scalar boson field theory) of the property of E under hermitian conjugation 



zE > (xi, x 2 ) 



zS > (x 1 ,x 2 ), (2.9) 



which holds for bosonic as well as (slightly modified) for fermionic field theories [22|. From 



both relations it follows that the difference E > (xi,x 2 ) — E < (xi,x 2 ) is real while the sum 
E > (xi,x 2 ) + E < (xi,x 2 ) is purely imaginary. 
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Furthermore, if the heat bath stays at thermal equilibrium at some fixed temperature 
and, as we have assumed, the self energy is solely determined by the property of the heat 
bath, the important relation 

£>(Jfe) = e fcf,/T £<(£), (2.10) 
holds, where E denotes the self energy in momentum space, 

E>«)(Jfc) = Jd\ Xl -x 2 )e ifc(xi - X2) £ >(<) (xi,x 2 ). (2.11) 
This relation follows e.g. from eqs. ( A19 ) when applying a Lehmann representation. Here 



we have implicitly assumed that the self energies stemming from the (isotropic) heat bath 
depend only on the difference of their space-time arguments. For the definition of the Fourier 
transformation the time integration is not limited by to but covers the whole range from 
— oo to +00. It is worth mentioning that our conventions are chosen such that iE < is always 
real and non-negative. In a transport theory it can be interpreted as the production rate 
for modes with the respective energy. 

For the following considerations it is useful to introduce in addition the retarded and 
advanced self energies 

E ret (x!,x 2 ) :=6(ti-f 2 ) [S > (x 1 ,x 2 ) -S < (x 1 ,x 2 )] = E++(xi,x 2 ) + E+-(x 1 ,x 2 ), (2.12) 
S av (x 1? x 2 ) := 0(t 2 - h) [S<(x!, x 2 ) - S > (x 1 , x 2 )\ = Z ++ (x u x 2 ) + Tr+( Xl , x 2 ) . (2.13) 



The property ( |2. 10| ) which is nothing but the Kubo-Martin-Schwinger boundary con- 



dition |£3[ will become crucial below when we will show that the modes of our system in 
the long-time limit will indeed equilibrate at the temperature of the heat bath. Explicit 
expressions for the self energies E, as given e.g. in appendix |A|, can be determined from a 
microscopic study of a system where some of its modes are assumed to be in thermal equi- 
librium and hence acting as a heat bath for the other modes. If these heat bath modes are 
integrated out one obtains an effective action for the remaining degrees of freedom, denoted 
with above. The influence of the heat bath modes is then contained in self energies and 
vertex corrections for the fields 0. Eq. ( p.lOj ) can be justified in the course of these calcula- 
tions. To keep things simple we restrict ourselves to self energies in the following and also 
neglect possible self interactions between the fields. 

As an example how the procedure described above works in detail we refer to ref. WA 



where the hard momentum modes (defined by some cutoff) of a scalar field theory are treated 
as a heat bath for the soft modes. Nevertheless we stress that the framework of the present 



paper is quite general. Instead of the specific example of ref. [H] one can also imagine 
e.g. a Yukawa theory where the fermions are integrated out or a system where all degrees of 
freedom except one act as a heat bath for the remaining one (cf. WW). 

Given the fields <p + and <p~ on the two branches of the CTP contour, four two-point 
functions can be defined. As usual we collect them as elements of a 2x2 matrix: 

D ab ( Xl , x 2 ) := -i{4> a {x x ) b (x 2 )> , a, b = +, - . (2.14) 

The expectation value in ( |2.14| ) is defined via path integrals as 
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(2.15) 



where p is the density matrix of the system for some initial time to and the normalization 
constant iV is given by 



N 



V\ 



e iS{<t>+,<t>-] p [ 



When evaluating the two-point functions by using ( |2.15|) one finds 

D ++ (x u x 2 ) = Q(ti - t 2 ) D>( Xl , x 2 ) + 6(t 2 - t x ) D < (x 1 ,x 2 ) 
D (xi, x 2 ) = Q(t 2 - h) D > (x 1 ,x 2 ) + Q(ti - t 2 ) D < {x x , x 2 ) 
D + ~(x 1 ,x 2 ) = -i(^(x 2 )0 + (xi)) =: D < (x 1 ,x 2 ), 
D~ + (xi,x 2 ) = -i((jT{xi) <p + {x 2 )) =: D > (x 1 ,x 2 ). 



The two-point functions are not independent of each other but satisfy plR 

D ++ + D— = D + - + D~+ . 
In addition, similar to eqs. (|2.8j), (|2.9|), the following relations hold 



D > (x 1 ,x 2 ) = D < (x 2 ,xi) , 

t < 
= iD > (xi,x 2 ) . 



iD > (xi,x 2 ) 

Again we now introduce retarded and advanced quantities by 

(xi,x 2 ) 



(2.16) 



(2.17) 
(2.18) 
(2.19) 
(2.20) 



(2.21) 

(2.22) 
(2.23) 



D TCt < 



D av (x u x 2 ) 



D ++ - d + - = ZT+ - d— 
e(h-t 2 ) [D>{x u x 2 ) -D<{x u x 2 )\ 
D ++ - D~ + = D + - - D~ 
e(t 2 -ti) [D<(x 1 ,x 2 )-D>(x 1 ,x 2 )} 



(2.24) 



(2.25) 



From ( p.l|) it is straightforward to obtain the equations of motion for the matrix elements 

D ab. 



(-□ -m 2 - 57 
(-□ -m 2 - 57 
(-□ - m 2 + 57 
(-□ - m 2 + 57 



)D ++ 
)D + - 
)D- 
)D- + 



£-+£>+- 

5T+D++ 



5, 
0, 
-5. 
0, 



(2.26) 
(2.27) 
(2.28) 
(2.29) 



where S is to be understood as a four dimensional 5-function with the appropriate argument. 
We note again that the "product" E.D is an abbreviation for / g? 4 x'E(;c, x')D(x', y). Since the 
four two-point functions are not independent of each other the four equations above contain 
redundant information. Therefore we will derive the equations of motion for retarded and 
advanced propagator and for an off-diagonal element, say 
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D < :=D + -. (2.30) 
Note that D ret and D av do not contain the phase space occupation number explicitly ET 



(see below). Hence the information about the time evolution of the latter is contained in 
the equation of motion for D K , only. Subtracting ( |2.27| ) from ( |2.26| ) and using ( |2.24| ) and 
( P-12j ) we can derive the equation of motion for the retarded propagator, 

(-□ - m 2 - £ ret )£> ret = 5 . (2.31) 

For the advanced propagator we get 

(-□ - m 2 - S av )D av = 5 (2.32) 

where we have used Q2.26Q , Q2.29Q , (|2.13|) and ( |2.7[) . One should remark, however, that 
the equation of motion ( |2.32| ) is completely equivalent to the one ( |2.31| ) for the retarded 
propagator and so does not contain any new information: From (|2.24 ) and ( |2.25|) and 
exploiting fl2.22| ) one has the familiar relation 

D ret ( Xl ,x 2 ) = D av (x 2 ,x 1 ) . (2.33) 

Additional information which is not already contained in ( |2.31D can be extracted e.g. from 
( |2.27| ) by expressing propagators and self energies in terms of quantities labeled with "ret" , 
"av" and "<": 

(-□ -m 2 )D < - £**£>< - S < J D av = 0, (2.34) 

which is nothing but the famous Kadanoff-Baym equation || . 

All equations of motion are explicitly causal as they contain only information from the 
past history of the evolution. The corresponding adjoint equations of motion of the above 
three equations follow directly by means of the relations (|2.8j), (pD|), ( p.22j ) and ( |2.23| ) 
and thus do not contain any further information. Hence, the equations Q2.31|) and (|2.34|) 
determine the complete and causal (non-equilibrium) evolution for the two-point functions. 

To get more insight in the (effective) action given in ( |2.1|) and the resulting equations of 
motion we introduce the following quantities: 

s(x 1 ,x 2 ) := ^sgn(ti - t 2 ) (U > (xi,x 2 ) - E < (x 1 ,x 2 )) , (2.35) 
a(xi, x 2 ) := - (S > (xi, x 2 ) - E < (xi,x 2 )) , (2.36) 
I(x 1 ,x 2 ) := -j (E > (x 1 , x 2 ) + E<(x 1 , x 2 )) , (2.37) 

which are all real valued due to the properties of S >!< discussed above. sgn(ti — t 2 ) denotes 
the sign function, i.e. 

1 for ti > t 2 

sgnih - t 2 ) := { for t x = t 2 } =20(^-^-1. (2.38) 



-1 for ti < t 



2 
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Recalling eq. ( |2.8| ) we find in addition that s and I are symmetric under the exchange of 
arguments, while a is antisymmetric: 



s(x 1 ,x 2 ) = s(x 2 ,x 1 ) , (2.39) 
a(xi, x 2 ) = —a(x 2 , xi) , (2.40) 
/((xi,a; 2 ) = /(x 2 ,a;i). (2.41) 

As we shall see in the following s yields a (dynamical) mass shift for the <fi modes caused 
by the interaction with the modes of the heat bath, a is responsible for the damping, i.e. 
dissipation of the fields while I characterizes the fluctuations. Especially we will show 
that the retarded and advanced propagator and hence also the spectral function of the (p 
modes are determined by s and a, only, while I only influences the number density. Note 
that the damping term is given by the difference of the transport theoretical gain and loss 
terms, S < and S > , respectively, while the fluctuation is characterized by the sum of both. 
This is quite plausible since there is no dissipation into other channels if gain and loss terms 
have the same size. On the other hand the fluctuations grow with any interaction. 

On account of ( gl35|) -(^37D together with we find that the CTP action (PD 

can be written as 



2 



-□ - m 2 ) <b + - 6 i-U-m 2 ) 



[8 + a) + d>-)+i 



(2.42) 



This expression is already identical to the influence functional given by Feynman and Vernon 
17 1 . (In appendix |A] we shortly review the original idea of how to obtain the influence 



functional when integrating out the bath degrees of freedom.) To the exponential factor in 
the path integral ( |2.15| ) the (s + a) term contributes a phase while the I term causes an 
exponential damping. Therefore we can indeed identify J as a 'noise' correlator. We will 
discuss that in more detail in section |TJ (see also appendix |A|). 

Now we can rewrite the definitions and the equations of motion given above in terms of 
s, a and I. Retarded and advanced self energies as defined in ( |2.12j ), ( ^.13| ) are given by 



S rct = s + a , (2.43) 

E av = s - a , (2.44) 

while is expressed as 

T 1 < = -a-il. (2.45) 
For the relevant equations of motion ( 2.31 )- (|2.34|) we find 



(-n-m 2 -s-a)D iet = d, (2.46) 
(- □ - m 2 - s + a)D av = 5 , (2.47) 
;-□ - m 2 - s - a)D < + (a + iI)D av = . (2.48) 
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We see that the last equation (the 'Kadanoff-Baym' equation) is the only one where / occurs. 
Hence the retarded and advanced propagators are determined by s and a while the number 
density is additionally influenced by the fluctuations /. 

The formal solution of the last equation including boundary conditions at the initial time 
t is 

oo 

D < (x l ,h; x 2 , t 2 ) = Jdt' dt" Jd 3 x' d 3 x" D™\x u t x ] x' , t') [-a - il}(x', t'; x", t") D aY {x", t"; x 2 , t 2 ) 



+ 



Jd 3 x' d 3 x" 



d 2 



D^ixuh; x 1 , t ) ——D<(x', f ; x", Q D™(x" } t ; x 2} t 2 ) 

- D™\x u tu x", t ) TrrD<(x", t ; x", t ) t '; x 2 , t 2 ) 

OIq OIq 

- ^r^ ret (^i, h\ t' Q ) W ,D<{x\ t ; x", t ') D av (x", t ; x 2 , t 2 ) 
+ Jr^ rCt (^i, *i! ^ Q D<(x', t ; x", t ) ^-jD av (x", t '; x 2) t 2 ) 

Ct OIq 



(2.49) 



-I io=* =*o 



We note that such a relation in the literature [^6|,[^,U] is sometimes denoted as a generalized 
fluctuation-dissipation theorem. We will see in a moment why. 

To get an intuitive interpretation of s and a we study the long-time behavior of this 
equation. In this case we can assume that the system becomes translational invariant in 
time and space and the boundary terms are no longer important. (Formally one might put 
to to — oo.) Then we get for the Fourier transform of D K the relation 



D<{k) = D Tet {k)[-a{k) - iI{k))D™{k) 



(2.50) 



while D and _D av are given by 



and 



D vet (k) 



D av {k) 



k 2 — m 2 — s{k) — a{k) 



1 



m? — s(k) + a(k) 



(2.51) 



(2.52) 



One customary way to study the influence of s and a is to look at the spectral function of 
the system. It is given by 



Aik) := - fl-(*)] = | fc2 _ m ^ 2 + |5W|2 

The spectral function fulfills the following energy weighted sum rule 

^ (k A{k, k ]^ 



(2.53) 



(2.54) 
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which follows as a direct consequence of the boundaries ( |2.24j ) imposed on the definition 



of D TCt (t = ti - t 2 < 0) = and the subsequent condition f t D TCt {t = t x - t 2 = 0+) = -1 
resulting from integrating the equation of motion ( |2.31| ) or ( |2.46| ) from t = 0~ to t = + . 
Note that a(k) is purely imaginary since a{x\,X2) is real and antisymmetric (cf. (|2.40|) ). 
For the same reason s is real. Inspecting the spectral function it becomes obvious that s 
contributes an (energy dependent) mass shift while a causes the damping of propagating 
modes, a is related to the commonly used damping rate T via 

f(Jfe) = 1 < ^L • (2.55) 

Let us now come back to the analysis of the equation of motion for D K . We introduce 
the long-time limit of the phase space occupation number n(k) via 

D<(k) = -2in(k)A(k) (2.56) 

and evaluate ( |2.50| ) using the spectral function: 

- 2i n{k)A{k) = D TCt {k)[-a{k) - iI{k)]D™{k) = - a ( fc )_~ t7 ( fc ) A (fy _ ( 2 .57) 

1 d[ Ki J 

To proceed we need a relation between a and I, i.e. a relation between the damping and 
the noise term. This is nothing but the fluctuation-dissipation theorem. In our case it is a 
simple consequence of the KMS condition (|2.10p using the definitions ( |2.36| ), ( |2.37| ): 



a(k)-il(k)_ £<(*) 1 g 



2a(k) ' E>(k)-E<(k) e k °/ T -l' 
Therefore we find 

— a(k) — il(k) 1 . . . 

" W = U) = s nM • (2 - 59) 

which indeed shows that the phase space occupation number of the soft modes in the long- 
time limit becomes a Bose distribution with the temperature of the heat bath. As we will 



briefly address in section [TV] the physical meaning of the KMS condition for the self energy 
is to guarantee detailed balance in kinetic collision processes in the long-time limit. 

If one now assumes that the coupling between the bath and the system becomes very 
weak, i.e. d,T,I — >• 0, the expression (|2.59| ) 

n(k) - ~ m ~ m (- ^ - 1 (2 60) 

n{k) ~ 2a(k) { '0' ) eVT _ i 

is still preserved as long as the KMS condition is fulfilled. (Starting from such a 'simple' 
observation we will comment in section on some similarities to the issue of the so called 
pinch singularity problem.) In the weak coupling limit the spectral function ( [2.53| ) becomes 
effectively onshell, 
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= ie-n*-$r + \mp ~ Isg,l(M Ke ~ m2) ' (2M) 



and D< thus reads 

- -27rzsgn(A;o)^(A: 2 -m 2 ) efco/T 1 _ i . (2.62) 

For our further discussion it is also interesting to explicitly write down the relation 
between a(k) and I(k): 

i ^= tU-^U im=eoai {^) iaik) - (2 - 63) 

In the high temperature (classical) limit we get 

I(k) = — 2i a(Jfc) , (2.64) 
k 

or, employing ( |2.55| ), 

I(k) = 2TT(k). (2.65) 

Recalling our discussion of Brownian motion in the introduction | this compares favorably 
well with (|1.8Q ! Indeed as we shall see in the next section we can define a quantity which 
obeys a Langevin equation very similar to the one for the Brownian particle. The physical 
meaning of I as a 'noise' correlator will become obvious. The relation ( |2.63|) thus already 
represents the generalized fluctuation- dissipation relation from a microscopic point of view 
by the various definitions of I, a and f through the parts E < and of the self energy. 

III. LANGEVIN EQUATION 

To see the connection between the formalism presented in section [TT] and stochastic 
equations we decompose the action S as given in ( |2.42| ) in its real and imaginary part and 
introduce the generating functional (see also appendix |X|) 



z[j+,r] ■.= jv[<i> + ,r]p[4>+A-}j s[ * + ^ ]+ij+ * ++ir *~ 

= Jv[(f) + , <f>~] p[(f) + , ^S[4>+ ) <t>~]+ij+4> + +ir4>--\ l ,4> + -4>-)i{4> + -4>-) 
TV J 

= 1 [v[Z]e-& J ~ l t fo[0+ <TM^ 01 e iRe5[ ^^ (3.1) 

/ V J J 



with 



N ■= [v^e-^ 1 1? . (3.2) 
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The generating functional Z[j + ,j~] in ( ft.ip can be interpreted as a new stochastic generat- 
ing functional Z'[j + + £, j~ — £] averaged over a random (noise) field £ which is Gaussian 
distributed with the width function J, i.e. 



z\j + ,r ] 



with 



and 



z'[j + ,r 



((0)):=lfvtOe-K'-H m 



(3.3) 



(3.4) 



(3.5) 



From the last definition we find that the (ensemble) average over the noise field vanishes, 
i.e. 



while the noise correlator is given by 



«a = o. 



(3.6) 



(3.7) 



which, in general, yields colored noise as long as / cannot be considered as being effectively 
white compared to the relevant time scales of the systems degrees of freedom. 

Obviously in the equations of motion generated by Z' the noise field £ acts as an addi- 
tional (randomly distributed) external source. Note that the action entering the definition 
of Z' is no longer S, but only its real part 



ReS 



(-□ - m 2 ) <p + - (-□ - m 2 ) 0~ + (0+ - 0") (-s - a) (<p + + <p~) 



-□ - 77T 



(-□ - m 2 



— (h a d) + d> a i 



(3.8) 



where we have used the symmetry properties ( 2.39 ) and ( 2.40 ) to derive the last expression. 
We introduce the following abbreviations: 



with 



(3.9) 



(3.10) 



J :-- 



(3.11) 
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and 



(3-12) 



-□ — m 2 — s —a 



R-=\ r, , 2 , • (3-13) 

Now we are aiming at a Langevin equation for a classical field. For this we start with 
the identity 

= Jv[<$>\ e i iM0+iJ *) = yb[$] p[$] + J) e ^ m+i J * . (3.14) 



Note that 

5 



?<5$(i,f 



= (3.15) 



for the case t > to we are interested in. Next we introduce the classical two component field 
where the average is defined as 



s j j z'[H 



(3.16) 



with an appropriate normalization N'. Obviously the equation of motion for ($)^ can be 
read off from ( |3.14 ): 



R($)t = -i. (3.17) 

This Langevin equation is indeed a generalization of the one for the Brownian particle 
discussed in the introduction. To see this we insert the definitions (|3.10| ) and ( |3.13|) and find 



(_□ _ m 2 _ s)(0±)? _ = _ e (3 lg) 

As at the initial time t the two fields and (<fi~)z do obey the same boundary conditions 
they stay the same for all times. This, of course, is what one would expect physically to be 
the case. Hence, for the expectation value 

<•)« = g) (3.19) 

holds. Finally, we thus get one equation of motion for <f)£-.=((j) + ) £=(<j)~) ^. 

(-□ -m 2 - s)0 e -a0 e = (3.20) 

Dropping all dependences on spatial coordinates this equation of motion becomes the ana- 
logue to ( |1 . 7| ) for the situation of a (quantum) Brownian oscillator instead of a free (and 
classical) Brownian particle: For a classical Brownian oscillator equation ( |1.7|) would be 
replaced by 
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f 

Mi + Mulx + 2 J dt'T(t-t')x{t') = . (3.21) 



(Similarly to the discussion in the introduction one can easily show that using relation ( |1.8|) 
in the long-time limit kinetic and potential energy become the same in accordance with the 
virial theorem and equipartition principle: ((p 2 ))/(2M) = Mu 2 ((x 2 })/2 = T/2.) With Q2.12 



and (|2.36| ) the spatial Fourier transform of the Langevin equation (|3.20|) takes the form 



<^(M) + (m 2 + k 2 )(f)^k,t) + J dt'Z rct (k,t-t')fa(k,t') = 

— oo 

t 

<j>t(k,t) + (m 2 + k 2 )(p^k,t) + 2 J dt' a(k,t - t') ^(k,t') = £(k,t) . (3.22) 

— oo 

Performing a partial integration we get 

t 

^(k,t) + (m 2 + k 2 -2T(k,At = 0))^(k,t) + 2 J dt' T(k,t - t') ^{k,t') = , (3.23) 

— oo 

Here we have defined the friction kernel T by means of 

J^r(M) := a(M), (3.24) 

and its Fourier transform in time does in fact correspond to the damping rate stated in 
equation (|2.55|) . In addition, in this form, the 'potential' part has been modified by a 
momentum dependent mass shift 

Am 2 = -2T(k,At = 0) = -2 [ —T(k,k ) = s(k,k = 0), (3.25) 

J 2n 



where the last relation follows from the definitions ( [2.55] ), ( |2.36|) and (|2.35| ). (From a purely 



mathematical point of view the transformation of the original equation of motion ( |3.2(J| ) to 
the form given in ( |3.23| ) can be interpreted as a subtraction method when expressing s and 
a of the retarded self energy by standard dispersion relations now with the involved friction 
kernel. We also note that in special applications the mass shift might be divergent and needs 
to be renormalized. This, however, is not important at this point of our considerations.) The 
analogy between the Langevin equation ( |3.23| ) and the one for the single (classical) oscillator 



( |3.21| ) is now obvious. The important difference, however, is the fact that the corresponding 
noise kernels / in ( |2.63j ) and ( |1.8| ), respectively, only agree in the high temperature limit. 
The interesting question now is to what extend the classical equations of motion (|3.18 



together with (K7h are an approximation for the full quantum problem given by the equation 



of motion (|2.48|) for D < . To figure that out we study the long-time limit of ( |3.17| ), i.e. when 
the boundary conditions at t can safely be ignored: 

($) f = _iTif. (3.26) 
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The inverse of R can be found via a Fourier transformation (in the long-time limit): 

k 2 — m 2 — s 



R 



—a 

—k 2 + m 2 + s 



R 



-i 



k 2 — m 2 — s —a 



[k 2 — m 2 — s] 2 + |a| 2 \ 
On account of fl3.12|) we find 



-k 2 + m 2 + s 



(3.27) 



(3.28) 



where _D ret is given by the Fourier transform of ( [2.51J) . Now we compare the full D K of ( 2.50 ) 
with 



The appearance of _D av is simply due to 

Jd 4 x' D Iet {x,x')Z{x') = Jd 4 x'Z{x') D av (x',x). 



(3.29) 



(3.30) 



Note that (|3.29|) is indeed the relation (|1.16| ) advocated in the introduction to hold in the 
(semi-) classical regime. In addition, using (|3.19p we realize that 



(3.31) 



Hence, taking equal time arguments, the expectation value of the following equal time 
commutator vanishes, i.e. 



*(([(0 + > € ,(<r> 5 l>>* 1= * 2 = o 



(3.32) 



which one would expect for a neglection of quantum corrections. More directly, comparing 
the classical approximation (|3.29|) with the exact solution (|2.50 ) we find that (—a — il) is 
approximated by —it. Of course this is justified, if 



\a\ < / 



(3.33) 



holds. Using the microscopic quantum version ( [2.63j ) of the fluctuation-dissipation theorem 
we get from (|3.33|) the condition 



ooth »1 



(3.34) 



Again we find that in the high temperature limit or - turning the argument around - for 
low frequency modes, 



k Q < T . 



(3.35) 



the classical solution yields a good approximation to the full quantum case. To be more 
precise: In simulations one has to solve the classical Langevin equation ( |3.17| ) and calculate 
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n-point functions by averaging over the random sources. This has been raised in [141. See 



also the work of |25fl for some practical examples regarding the mathematical evaluation of 
expectation values in the transition from quantum to classical field theories using dimensional 
reduction techniques. 

We now go one step further and study the equations of motion for the quantum two-point 
functions with external noise. We introduce the 'noisy' propagators 



Df(x 1} x 2 ) :=-t((t> a {x 1 ) ( t>\x 2 )) i 



(3.36) 



and derive their equations of motion by differentiating ( J3.14Q with respect to J and setting 
J = £ afterwards: 



(3.37) 



Since the 'noisy' propagators are extracted from a closed time path action (with ReS instead 
of S) they satisfy the relation (cf. |2.21| ) 



Again we introduce the quantities 



d; 



d: 



(3.38) 



Df 
Df 



DP 
Dp 



D: 



From the components of ( 3.37 ) we getQ 



m 



D 

D: 



a)Df = 5 , 



(3.39) 
(3.40) 
(3.41) 



(3.42) 



m 



a)D 



av 



and 



(3.43) 



-□ 



m 



- s-a)Df + aD? = i£. 



(3.44) 



They ought to be compared with ( 2.46| )-( p.48 ). We immediately recover that and -D| v 
obey the same equations of motion as _D ret and D av , respectively. Since also their respective 
retarded or advanced boundary conditions are the same we conclude that 



Df 
Df 



D mt 
D av . 



(3.45) 
(3.46) 



Only the equation of motion for D^ and hence for the occupation number is modified 
compared to the one for D<. One could have already guessed that from the equations of 



L We set (cj) + )^=((j) }£=0£, cf. discussion after eq. ( p. 18 ). 
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motion (|2.46|) -( p.48|) since only in the last one of these equations the noise kernel I appears. 
The latter, however, makes the difference between the generating functionals Z and Z' 
(defined in ( |3.1| ) and ( |3.4| ), respectively) from which the respective equations of motion are 
deduced. 

To rederive the Kadanoff-Baym equation ( |2.48| ) we have to average ( |3.44| ) over the noise 
fields. To calculate ((£ <p^j) we solve the equation of motion ( |3.20| ) for (p^. 



to 



(3.47) 



+ d 6 x 



3 ' 



-Z7 ct (x, t; x' , t' ) Jr<M^ t' ) + Jr^ rot (x, t; x', t' ) tjt&t , to) 



t =t' 



Using ( |3.6| ) the terms caused by the boundary conditions drop out from the calculation of 
((£ an d we find 



-ID 3 



(3.48) 



where we have used (|3.30|) to change D ret into D av . Thus averaging (|3.44|) over the noise 
fields according to (|3.5| ) we indeed rederive the Kadanoff-Baym equation fl2.48|) . Hence, as 
mathematically expected by inspection of ( |3.3j ), (ensemble) averaging the 'noisy' propagator 
D? b over the quantum noise, just yields 



({Df(x 1 ,x 2 ))) = D ab ( Xl ,x 2 ). 



(3.49) 



This shows that the Kadanoff-Baym equation can be interpreted as an ensemble average 
over fluctuating fields which are subject to noise, the latter being correlated by the sum of 
self energies S < and S > , i.e. from a transport theoretical point of view the sum of production 
and annihilation rate. We will discuss the meaning of eq. ( |3.44|) in more detail in the next 
section. We close this section by noting once more that the 'noisy' or fluctuating part 
denoted by I inherent to the structure of the Kadanoff-Baym equation ( 2.34Q guarantees 
that the modes or particles become correctly (thermally) populated, as can be realized by 
inspecting ( |2.59| ) or ( |2.60| ). / is related to the usual damping rate ( |2.55| ) by means of a 
generalized quantum version (|2.63|) of the fluctuation-dissipation relation. 



IV. TOWARDS A MICROSCOPIC KINETIC BOLTZMANN-LANGEVIN 

EQUATION? 

In this section we want to elucidate more on the physical content of the Kadanoff- 
Baym equation(s) ( |2.31| ) and ( |2.34j ) and the corresponding equation ( 3.44 ) for the 'noisy' 



counterpart. One major goal is to derive a standard kinetic transport equation for the (semi- 
classical) phase-space distribution f(x,k,t) of bosonic quasi-particles from the equation of 
motion (|2.34|) which should be valid in a weak coupling and nearly homogeneous system. 
Various derivations are already presented in the literature [p|p|, p2| , p6| , p7| , so that we will 
put some more emphasis on the validity of the necessary steps and assumptions which one 
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typically has to make. The resulting kinetic transport equations have to be understood as an 
ensemble average over possible realizations of the system. On the other hand, one can also 
try to pursue a completely analogous expansion for the noisy counterpart, eq. ( |3.44| ). The 



resulting equation resembles in its structure the phenomenologically inspired Boltzmann- 
Langevin equation, which has the form of a standard transport equation supplemented by 
an additional noisy (i.e. fluctuating) term. This term allows for a quantitative description 
also for the (thermal) fluctuations of the phase space distribution. A closer inspection of our 
equation, however, shows that its interpretation does not exactly correspond to that of the 
Boltzmann-Langevin equation. It actually will not only describe the thermal fluctuations 
of the phase space distribution. In fact, the suggested interpretation of the introduced 
function f^(x,k,t) as an individual phase-space distribution is not correct, although its 
ensemble average does yield the average phase-space distribution. In turn, the validity of 
the phenomenological Boltzmann-Langevin approach may be doubted. 

Before we explicitly address the issue of a kinetic transport description we first will study 
the average particle number (n(k,t)) and its fluctuation (n 2 ) — (n) 2 at some given late time 
t, when the system stays at equilibrium. This will help to understand the later procedure 
and its interpretation. 

The (standard) second quantized description for an (interacting) real boson field at some 
given time t is given by 

r d^k 1 

<f>(x,t) = J — -j= (a n (t)e~^e ik * + a^e^e" 4 ") . (4.1) 

a^(t) denotes the annihilation and at(t) the creation operator of a boson of momentum k 
at time t. It is assumed that the quantized modes at that particular time t are propagating 
to some major part with some (appropriately chosen) frequency ko = u^. This corresponds 
to a quasi-free propagation of the bosonic particles. As the particles are interacting, i.e. in 
our special case with the environmental heat bath, the creation and annihilation operators 
still do possess an explicit time dependence. The separated, oscillating phases e^^ 1 thus 
correspond to the representation of the fields within the interaction picture. For the field 
momentum one has 

7r(£,t) = J^(£,t) (4.2) 

J-L j^l (_ ia (+y-^tj%S + i Utyi^t -ikl\ 

(27r) 3 V 2 V fcW fcW > 

In obtaining the second line the explicit time dependence of the creation and annihilation 
operator has been neglected. This approximation can thus only be valid if the dominant 
time evolution of the field operator is described by the phase e £*. Physically speaking, the 
validity of this approximation corresponds to the validity of the quasi-particle approximation, 
i.e. to a quasi-free propagation of the bosonic particles. In the following discussion we restrict 
ourselves to this approximation. 

Solving ( |4.1| ) and (|4.2|) for and at, one has 
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Tf0(-M) 



1 -9 Ll r x 
t—(p{—k, t) 



2ut dt 



(4.3) 



For the particle number (N(k,t)) of particles with the momentum k inside the system at 
the given time t one now writes 



(AT(M)> = <4 a *> 



(4.4) 



u h 1 d d i . d d A . , . -* . , - .. . 

— + 7- — — — + -(- — ) (6(-k,t)6(k,t')) 

2 2u e dtdtf 2 K dt' dt ! r ' m ' ;/ 



t'=t 



The two-point expectation value can readily be expressed by means of the off-diagonal entry 
D < of the CTP one-particle propagator as 



(0(-M)0(M')> 



t,t -+oo 



V 



d 3 x / d 3 y e^e"^ {iD< (y, t'; x, t)) 



(iD<(k, ko) 



dk 



£ Jk (t-t') 



2tt 



(4.5) 



where the second line only holds in the long-time limit when the system is at equilibrium 
(and assumed to be homogeneous). Combining (|4.4|) and ( |4.5| ) with (|2.56 ), the particle 
number density (n(k)) = (N(k))/V becomes 



r dk ( k 2 \ 
(n(k)) = J I u % + ^ + 2k j A(k, k )n B (k ) . 



(4.6) 



In the quasi-particle regime where the approximation carried out in ( |4.2j ) is valid one assumes 
for the (nearly onshell) width 



r(fc, k « u%) < wjj = \/k 2 + m 2 , 
so that the spectral function becomes quasi onshell (see (|2.61 )) 



A{k) « 7rsgn(A;o) 5(k 2 - m 2 
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k 



(4.7) 



In this limit the appropriate choice for is then obviously given by uj^ = One thus 
ends up with the expected result 



(n(k)) « n B {4) , 



(4.8) 



i.e. the average particle number density is given by the thermal Bose distribution evaluated 
at the quasi-particle energy u^. 

For our further discussion it is illustrative to write down the corresponding particle 
number (N(k, t))% resulting from the 'noisy' generating functional Z'[j + + £, j~ — £] of ( 13.41) , 
or, more directly, from the 'noisy' off-diagonal propagator D^: 



(N(k,t))t = (oW> f 



(4.9) 



UJt 



1 d d i , d d 



2 ' 2u n dtdt' 2 y dt' dt' 



(0(-M)0(M')>< 



t'=t 
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with 



(4.10) 



One is tempted to interpret the such calculated number (N(k,t))^ as the particle number 
in the momentum state k for one particular noise realization contributing to the overall 
generating functional by means of ( |3.5| ). In the long-time limit, is formally given by 
means of ( |3.44|) and (|3.47| ) 



-i(D ret £)(£D av ) - D mt aD 



(4.11) 



Due to the inherent fluctuations the propagator will depend explicitly on each coordinate 

separately. Because of the property (|3.49|) the ensemble average of (N(k,t))^ is, of course, 
just the average particle number ( [4.4|) , i.e. 



(«iV(MM = (N(k)) . 



(4.12) 



Hence, (N(k,t))^ will fluctuate around its average value. One would thus be tempted to 
think of its physical interpretation as the fluctuating particle number. However, as we will 
see below, this will not be the case. 

The average fluctuation in particle number is defined via 



(N 2 (k,t)) - (N(k,t)) 2 = (a^a^) - (ata u ) 



The first expectation value can now be obtained in a similar fashion as 
making use of (43) one finds 



(4.13) 
Repeatedly 



( a l a k a l a k) 



1 

4 

+ 



lot + 



d d 



d_ d_ 

dti dt 2 



d_ d_ 

dt 3 dt 4 



d d d d d d 
dti dt 2 dt 1 dt 3 dti dt 4 dt 2 dt 3 



d d d d 
dt 2 dt 4 dt 3 dt 4 



% 

+— 



d d d d d d d d d d d d 
dt\ dt 2 dt 3 dti dt 2 dt 4 dt 2 dt 3 dt 4 dti dt 3 dt 4 



+ 



d d d d 



uji \ dti dt 2 dt 3 dt 4 



((/)(— k, ti)(p(k, t 2 )4>(-k, t 3 )(j)(k, 4 )) 



t\=t2=t3=t4=t 



(4.14) 



In our special case the effective action ( |2.1| ) is only quadratic in the fields. Therefore the 
four-point function can straightforwardly be evaluated by means of the generating functional 
(^.1[). It reads (with a, b,c,d = +, — ) 



{<t> a {xi)<l> h {x2)<t> C {xz)<t> d {^)) 



(4.15) 



(x 3 )</> d {x 4 ) e lS ^ + ^~ 



D ab (xi,x 2 )D cd (x 3 ,x 4 ) + D ac (xi,x 3 )D bd (x 2} x 4 ) + D ad (x u x 4 )D bc (x 2} x 3 ) 
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This factorization corresponds to the usual Wick decomposition valid for Gaussian function- 
al. For the correct ordering one considers the limit ^1,^2,^3,^4 — * t with t\ > t 2 > t 3 > £ 4 : 

(<f>{-k, h)<t>{k, h)<f>{-k, t 3 )(f>{k, U)) = (4> + {-k, t 2 )<p + (-k, t 3 )<p + (k, U)) . (4.16) 

By using ( |4. 15| ) and rewriting the appearing contractions as (using ( |2.17| ) and ( |2.22| )) 

D ++ (xi, Xj) = D > (xi,Xj) = D < (xj,Xi) for U>tj , (4-17) 

the four-point function entering into ( |4.14j ) can be expressed by a sum of products of two 
D < propagators. Now one can explicitly evaluate ft4.14|) , which is straightforward but cum- 
bersome. 

Instead, one can also decompose the expectation value {^\a^a}^a^j directly by employing 
the (just derived) Wick contraction property as 

( a i a k a i a k) = ( a l a k)( a l a k) + ( a l a k)( a k a l^ + ( a l a l^ a k a k) ■ ( 4 - 18 ) 

One can easily show that (atop) = (c^g) = 0. For (a^a^) one finds a similar expression to 
([4.4|)-(p~6|) for the long-time limit 

r dkn ( k 2 \ - 

(«a?4) = I 2^ [ u k + ^ - 2k °) Ak, k )n B (k ) . (4.19) 

At this stage we note that in general for a strongly interacting theory where the damping is 
large we find 

^<a«4> ^ 1 + (420) 

for a spectral function A differing from a quasi onshell form given in ( |4.7|) . This basically 
means that the particle number cannot be unambiguously defined for a spectral function 
which is not (nearly) onshell. 

Only in the quasi-particle limit (cf. ( [4.71) and ( |4.8|) ) which we consider here, ( |4.19| ) be- 
comes 

i(a fc -4) « 1 + n B (ul) . (4.21) 

Putting everything together, the familiar result for the average (thermal) fluctuation is then 
recovered: 

(n\k)) - (n(k)) 2 « n B (wJ) (l + n B (cJ§) . (4.22) 

It is valid in the quasi-particle regime at equilibrium. 

If the interpretation of (N(k,t))^ as a fluctuating particle number would be correct, we 
should also get the rhs. of (|4.22|) by calculating the quantity 



^{(((N(lt))l)}-(((Nfy))^) . (4.23) 
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In the following we will first derive a transport equation for the noisy occupation number 

(n(k,t))t:=±{N$,t))t. (4.24) 

We will indeed find an equation which resembles a Boltzmann-Langevin equation support- 
ing again the idea that (N(k,t))^ might be interpreted as a fluctuating particle number. 
Afterwards we will check whether ( |4.23| ) is indeed identical to ub^b + !)• 

To derive kinetic transport equations for the occupation number as well as for its noisy 
counterpart on the same footing we first discuss the derivation of the former. We start 
with the exact equations of motion ( p.46j )- p.48| ) for D iet , .D av , and D < . If the system is in 



a general off-equilibrium state these two-point functions depend not only on the so called 
microscopic variable u := X\ — x 2 but additionally on the macroscopic center-of-mass variable 
X := (xi + Xz)/2. In a transport approach one generally assumes that the dependence of 
the propagators on the latter is rather weak. Note that for the equilibrium case the two- 
point functions do not depend on X at all. In any case the propagators are rapidly varying 
functions with respect to the microscopic variable u and it is therefore useful to perform 
a Fourier mode decomposition with respect to u. This decomposition is called the Wigner 
transformation. For our case at hand we have e.g. for D < \ 

D < (X,k) = Jd\ie iku D < (X + u/2,X -u/2) . (4.25) 

(Note that here the time integral is not restricted by t . One typically argues that one puts 
to —>■ — 00 in order to define the Wigner transformed quantities. A finite to would result in 
some switching-on artifacts.) To transform the equations of motion ( |2.46| )-( f2.48D we have to 
calculate the Wigner transform of convolution integrals. This is presented in appendix ||. 
The result for a convolution of two so far arbitrary functions £ and D is 

jdSi e lku JdSj E(X + u/2, y) D{y, X - u/2) = e - i0 £(X, k) D(X, k) (4.26) 



with the operator O being defined as [£8 

1 



2 



d^-d^ . (4.27) 



We start with the equations for retarded and advanced Green's functions. Applying the 
transformation property ( f4.26| ) to the convolution of (—O xl —m 2 ) 5(xi—y) — s(xi, y)^fa(xi, y) 



with D Ict / SLV (y,x 2 ) we find 

e- t0 [k 2 - m 2 - s{X, k) =F a(X, k)]D Ict/! " v (X, k) = 1 . (4.28) 

In principle, in our case the self energies s, a, and / do not depend on the macroscopic 
variable X since their Fourier transforms depend on the difference variable, only. The 
reason for this is our basic assumption that the self energies are generated by an isotropic 
heat bath (cf. ( [2.11| ) and the remarks thereafter). In general, of course, also the self energies 



might depend on the macroscopic variable X. This holds generally for the more standard case 
where the self energies are generated by the self interaction of the modes under consideration. 
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In spite of the fact that we are dealing with a simpler case here we want to stress that the 
following considerations remain valid in the more general case where the self energies do 
depend on X. Therefore we keep the dependence on X also for the self energy functions. 
Next we separate (|4.28|) in real and imaginary parts. To do so we first note that 



(D 1 



D B 



ReD ret + iA 



(4.29) 



which can be easily shown by using the definitions for retarded and advanced Green's func- 
tions (|2.24|) and ( |2.25|) together with relation ( [2.22|) . As already defined in ( |2.53| ) for the 
special case of a homogeneous (i.e. X independent) system, the spectral function in its 
Wigner representation is given by 



A(X,k) :=-[D™\X,k)-D™(X,k)\ 



lmD av (X, k) = -ImD Te \X, k) . 



(4.30) 



In addition we recall that s is real while a is purely imaginary which can be seen by Fourier 
transforming the symmetry relations ( |2.39| ), ( |2.40| ). Inserting ( f4 . 2 9| ) in ( [4.28| ) we get the 
following set of equations 



cos O 


[{e 


— m 2 


sin O 




— m 2 


cos O 




— m 2 


sin O 


[{*> 


— m 2 



- s} {ReD TCt } + {ia} {A} 

- s} {lW ct } + {ia} {A} 



1, 
0, 
0, 




(4.31) 
(4.32) 
(4.33) 
(4.34) 



where the pair of functions the operator O acts on are now put in curly brackets. 

So far the equations of motion are still exact. Now we make use of the basic assumption 
of transport theory that all quantities depend only weakly on X. Therefore higher orders of 
the operator can be neglected. The expansion in orders of is called gradient expansion. 
Transport equations are derived by considering terms up to linear order in 0. Within this 
approximation scheme we find 



(k 2 



m 



s)ReD Ict + iaA 







{k 2 - m 2 - s} {ReD ret } + {ia} {A} 
(k 2 - m 2 - s)A - iaReD TI 
[k 2 -m 2 -s) {A} - {ia} {ReD rct } 



1 






(4.35) 
(4.36) 
(4.37) 
(4.38) 



The purely algebraic equations ( |4.35|) and (|4.37|) can be readily solved yielding results well- 
known from the homogeneous case: 



A{X, k) 



i a(X, k) 



[k*-m 2 -s(X, k)} 2 + \a(X, fc)p 



and 



ReD TCt (X, k) 



k 2 



m 



s(X, k) 



[k 2 -m 2 -s(X } k)] 2 + \a(X, k)\' 



(4.39) 



(4.40) 
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or in terms of retarded and advanced Green's functions (cf. ( |2.51| ), ( [2.52| )) 



D ^ x ' k) = e- m *- t( x, k ) + nx,k) ■ (442) 

An important thing to note here is that these relations hold not only in zeroth order of the 
gradient expansion (homogeneous case) but are also still valid in linear order. Actually the 
relations ( |4.36| ) and (|4.38|) provide no additional information. A tedious but straightforward 
calculation shows that A and Re.D ret as given in ([4.39|) and ( |4.40| ) identically fulfill fll36|) 
and flPBD . 

Having especially calculated _D av up to linear order in O we can proceed and study the 
transport equation resulting from ( |2.48| ). After Wigner transformation we get 

e- l °[k 2 - m 2 - s(X, k) - a(X, k)]D < {X, k) = - e -' J °[a(X, k) + iI(X, k)]D av (X, k) . (4.43) 

We note that D < is purely imaginary which can be checked using relation ( |2.23| ). Guided 
by the equilibrium case ( |2.56|) we introduce the (real valued) occupation number n(X, k) by 
the following ansatz (also cf. f27| , |2~8|l ): 

D<(X, k) = -2in(X, k)A(X, k) . (4.44) 

Similarly we decompose a + il by introducing the self energy occupation number n^: 

a(X, k) + iI(X, k) (= -£< (X, k)) = -2a(X, k)n^ (X, k) . (4.45) 



For the special self energies we have introduced in section 0, the self energy occupation 
number is simply given by the Bose distribution as already calculated in ( |2.58| ). However, 
also for arbitrary self energies it is possible (and, as we shall see below, useful) to introduce 
ns. Since a is purely imaginary while / is real we find that uy, is also real. Decomposing 
( |4.43| ) in real and imaginary part and keeping terms up to linear order in the gradients yields 

ianA - 0{k 2 - m 2 - s}{nA} = n^iaA - 0{ian s }{ReD iet } (4.46) 

and 

(k 2 - m 2 - s)nA + 0{ia]{nA} = ian^ReD Ict + 0{ian^}{A} . (4.47) 

From the definition of the gradient operator O in fl4.27p we find the useful property 

0{A}{BC} = BO{A}{C} + 0{A}{B}C (4.48) 

for arbitrary functions A, B, and C. (Actually O is a four-dimensional generalization of 
the Poisson bracket well-known from classical mechanics.) Using ( |4.38| )-( fl.40| ) and (|4.48|) we 
find a transport equation 

ianA - 0{k 2 - m 2 - s}{n}A + nO{ReD vet }{ia} 

= tan^A + 0{ReD iet }{nx}ia + n E 0{ReD ret }{ia} (4.49) 
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and a generalized mass shell equation 



(k 2 - m 2 - s)nA + 0{ia}{n}A + nO{ia}{A} 

= (k 2 -m 2 - s)nxA - 0{A}{nx}ia + n^O{ia}{A} . (4.50) 

Readers who are familiar with transport equations e.g. of Vlasov or Boltzmann type might 
feel uneasy with the unillustrative form of Q4.49 ). At first sight it might even be unclear 



where the obligatory drift term is hidden. At least that we will show right now before 
introducing additional assumptions under which more familiar expressions for the transport 
equation can be derived. The drift term is generated from the gradient operator acting on 
the k 2 term: 

0{k 2 }{n} = -k^n. (4.51) 

Without additional assumptions the transport and generalized mass shell equation cannot 
be simplified any further. If we assume, however, that the damping of the modes is small 
it should be possible to recover the standard quasi-particle picture of (kinetic) transport 
theory. To incorporate that we drop all terms in ( |4.49| ) and ( 4.50|) which contain both the 



gradient operator and the damping a. Note that in the spirit of the gradient expansion the 
terms linear in O are already small compared to the o(O ) contributions. If additionally 
a is assumed to be small it should be justified to neglect terms ~ aO. We also neglect 
terms where s appears together with O since s and a are generated from the same self 
energies (cf. ( |2.35| ), ( |2.36| )) and are connected by a dispersion relation. With the additional 
assumption of a being small we get the much simpler equations 

- 0{k 2 - m 2 }{n}A = ia(nx - n)A (4.52) 

and 

(k 2 - m 2 - s) (n - n E )A = . (4.53) 



For small a the spectral function ( f4.39| ) becomes strongly peaked at the mass shell. This 
corresponds to the quasi-particle regime. Hence, then (k 2 — m 2 — s)A (approximately) 
vanishes and the last equation is automatically fulfilled no matter how n and are related. 
This motivates the ansatz (|4.44|) . In addition the appearance of the spectral function in 
the transport equation ( |4.52j ) effectively puts n and ns on the mass shell. To rederive the 
commonly known collision term we only have to insert the definitions for and ia: 

ia(n s -n) = ^(E> - E<) ~ n ) = I S< ( n + 1 )~ ^ >n ) ■ ( 454 ) 

( [4.6D suggests to define now the onshell phase-space distribution function f(x,k,t) in the 
quasi-particle regime as 

-^Ak A(X,k,k )n(X,k,k ) « n(X,k,k ) Q . (4.55) 
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As in ( fO|) the factor 4ko takes into account the normalization of the spectral function stated 
in ( 2.54Q . The volume average of the so defined distribution function just gives the average 



particle density of momentum state k, i.e. 

i J d 3 xf(x,k,t) = (n(k,t)) , (4.56) 

whereas the integral over momentum states results in the local particle number density 

d 3 k 



(2n) 



f(x,k,t) = (n(x,t)) . (4.57) 



oo 

dko, 



Integrating the equation (|4.52|) with / -^-^ko . . . and using ( [4.54| ) we find in the end 

the well-known relativistic kinetic transport equation in the quasi-particle approximation 
(cf. e.g. ©0) 



k^f(x,k,t) = l Ux < (X,k)[f(xXt) + l]-iX>(X,k)f{x,k,t)) h _ , (4.58) 

K 

where all energies have to be evaluated onshell. 

While this is the commonly used form of the transport equation the unfamiliar form 
( |4.52j ) allows for an interesting interpretation (also cf. ||28|| ). Obviously it has the form of 
a relaxation equation where the relaxation time is given by the damping rate T as can be 
easily seen by concentrating on a system which is homogeneous in space. In this case the 
transport equation ( |4.58|) reduces to 



^/(M) = -r(i,fc, W ?)(/-/ E ) , (4.59) 



where /s is defined analogously to (|4.55|) and we have used the (generalized) relation (|2.55|) 



fc r(X, k) = ia(X, k) . (4.60) 

Furthermore we note that the collision term is proportional to the difference between the 
actual occupation probability / and the one demanded by the self energies, f%. This term 
drives the system towards equilibrium and only there the two quantities / and /g become 
identical, i.e. they both approach the onshell Bose distribution. 

For the self energies we have introduced in section || which satisfy the KMS condition we 
find that is already identical to the Bose distribution. In this case ( |4.59j ) is identical to the 



traditional relaxation equation which is usually set up if small deviations of the occupation 
number from its equilibrium value are studied. 

Finally, we would like to mention briefly how the simplest kind of transport equation 
emerges in the formalism presented here. If we had included a Hartree contribution in our 
self energies the mass term m 2 would be replaced by a space-time dependent effective mass. 
Neglecting all other self energies (S > ' < ) we find the Vlasov equation 

- 0{k 2 - [m(X)] 2 }{n(X, k)} = k„d^n(X, k) + d%m(X)d*n(X, k) = . (4.61) 
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Naively one might have expected that the Vlasov term is generated by the mass shift s. 
This however cannot be true since according to ( 2.35Q the latter is generated from the self 



energies S >,< which do not include the space-time local Hartree term (cf. e.g. p2| , |2"9fl ). 

From our rather lengthy derivation it should become clear that the use of a kinetic 
transport description is no longer an adequate description, if 

• the systems degrees of freedom do not lie in the quasi-particle regime (the situation 
encountered in a strong coupling theory), i.e. the spectral function A cannot be ap- 
proximated by a nearly onshell form, 

• the system reacts violently in the sense that a clear separation between the intrinsic 
scale u and the macroscopic scale X does not really exist, i.e. a standard Markovian 
approximation is not really justified. 

One could easily argue that such a situation is met in describing (relativistic) heavy ion 
collisions. However, due to the complexity, more or less all numerical descriptions of the 
transport processes are based on the semi-classical kinetic limit. For first numerical inves- 
tigations trying to go beyond that limit for a fermionic system (a system of interacting 
nucleons) we refer the reader to p0|-p2|. 

The kinetic transport equation ( |4.58| ) was obtained by a systematic first order gradient 
expansion of the underlying Kadanoff-Baym equations ( |2.46| )-( ]2T4"5| ). It has to be understood 
as the evolution of the ensemble averaged quasi-particle phase-space distribution /. As we 
have stressed at the beginning of this section, this distribution will of course fluctuate around 
its ensemble average. Can one also describe the evolution of the phase space distribution 
including the inherent fluctuations? One is immediately tempted to redo our above analysis 
for the noisy counterparts ( |3.42| )-( |3.44D of the Kadanoff-Baym equations derived in section 



IJ. Due to the identities (|3.45|) and (|3.46|) the spectral information for the noisy counterpart 



is exactly the same. In turn the same expressions ( |4.39| ) and fl4.40p will emerge and do not 



depend on a particular noise sequence. It is sufficient to concentrate on the noisy transport 
equation ( |3.44| ). For this we rewrite Q3.44 ) by employing ( |3.47| ). In the latter equation we 



neglect the (unimportant) second part which takes care of the explicit initial conditions for 
the fluctuating field being specified at some early time t . One then has 

(-□ - m 2 - s - a)D< + (a + il) Df = - I) Df := -i(AI) Df . (4.62) 



In this form the lhs. is completely equivalent to that of (|2.48| ). The rhs., however, is now 
substituted by a fluctuating source 

A/ , «AI» = (4.63) 

which vanishes on average. 

We now again assume that all quantities are smooth in their macroscopic center-of- 
mass variables X, especially also and AI. One should note here that this assumption 
is not obvious. For an arbitrary noise sequence there is a priori no reason why £(X + 
u/2)£(X — u/2) should depend only weakly on X. Since the noise terms influence the 
quantity D^(X + u/2,X — u/2) the same remark holds true for the latter. On the other 
hand, for a transport theoretical approach to be reasonable the noise correlator I(X + 
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u/2,X — u/2) has to depend only weakly on X like all other self energies. (Otherwise the 
gradient expansion performed above would be invalid.) Therefore one might assume that 
the product of noise functions where the dependence on X is strong do not dominate the 
stochastic process under consideration. Without that assumption it would be impossible to 
derive a transport equation for the noisy counterpart of the phase space density since the 
gradient expansion could not be performed. 
The Wigner transform of AI is defined as 



AI(X, k) = JdSi e iku (£(X + u/2)i(X - u/2) - I(X + u/2, X - u/2)) 



AI(X, fc) y 
Similar to fl4.44|) we make the ansatz 

Df(X, k) = -2inz(X, k)A(X, k) 



(4.64) 



(4.65) 



The further analysis is now completely analogous to the one given above. The corresponding 
equations to (|4.52| ) and ( 4.53|) take the form 



and 



- 0{k 2 - m 2 }{n^}A = ia(n s - n^A + ^-A 



(k 2 -m 2 -s) \nt - n s - — ) A = 
\ 2ia I 



(4.66) 



(4.67) 



By the arguments given after ( |4.53|) this last equation again does not contain any new infor- 
mation as (within the quasi-particle approximation) it is automatically fulfilled. Finally, we 
define the following quantity f^(x, k,t) (which one would like to interpret as the fluctuating 
onshell phase-space distribution) as 



iik 

ft(x,k,t) := / -^4k A(X,k,k )n^(X,k,k ). 



(4.68) 



The kinetic transport equation incorporating the inherent fluctuations then results as 



k^fcix, k,t) = \ (it<(X, k) k, t) + l}- iY>(X, k) ft(£, k, t)) 



k =uj° 

k 



(4.69) 



+ ujI^(x, k,t) 



where all energies have to be evaluated onshell, and the fluctuating source is introduced 
via 



^(Xjkjt) :- 



k o 



dko 
~2^ 



2k A(X, k, k )AI(X, k, k ) . 



(4.70) 



Dropping all the x- and fc-dependences, ( |4.69| ) has the more intuitive form 
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j t m = -f (m - h) + w) • (4-7i) 



The just derived kinetic transport process Q4.69Q and ( 4.71Q has the structure of the 



phenomenologically inspired Boltzmann-Langevin equation [Q. In order also to describe 
fluctuations around the average, Bixon and Zwanzig [[33[ postulated that in analogy with 
the Langevin equation for a Brownian particle, these fluctuations should be described by 
a stochastic (classical) Boltzmann equation. Such a stochastic kinetic equation has indeed 
been obtained by adding a fluctuating collision term to the linearized Boltzmann equation 
(of the form of equation ( |4.58| )), i.e. it is of the form 



j t h{t) = -rW) - /„) + rf ea (t) , (4.72) 



for a system near equilibrium. The correlation function of the fluctuating collision term 
was derived on the basis of the fluctuation-dissipation theorem. Also it was shown that 
the stochastic Boltzmann equation provides a basis for describing hydrodynamic fluctua- 
tions. On the other hand, the equations of Bixon and Zwanzig were not derived from first 
principles but obtained on the basis of intuitive arguments. In the following we translate 
their arguments to our framework and try to motivate how the correlation function of the 
fluctuating source in (|4.72|) accordingly should look like. Afterwards we will compare this 



intuitively 'derived' results with the exact microscopic calculation. 

In its form the equation (|4.72| ) resembles the form of the Langevin equation ( |1.1|) in 



the introduction, especially if one introduces the substitution fa = f eq + A/| and writes 
down the resulting equation for the fluctuating quantity A/g. At this point, however, the 
correlation ((JF| hen (t).7 : |' hen (£'))) is not yet specified. One now simply postulates that jF| hen 
is described by a Markovian (i.e. white) and Gaussian stochastic process of the kind similar 
to ( |1.2|) , i.e. one takes 

((J-f en (t))) = (4.73) 
((jf^Wjf^tf))) = lf cn 5{t - t') . (4.74) 

In the long-time limit, similar to ( |1.4j ), the average in the square of A/t is then readily 
obtained as 

rphen 

(((A/,) 2 )) ^ ^r- • (4-75) 

In our case, for Bose (quasi-)particles, and recalling the familiar result ( (4.22[ ), (((A/e) 2 )) 
should describe the thermal fluctuation in the particle number density, i.e. (with / eq = n B ) 

(((A/ 5 ) 2 » := n B (l+n B ) . (4.76) 

Hence, the correlation strength J^- hen of the fluctuating source, is suggested to be given by 

lf cn = 2f n B (l+n B ) . (4.77) 

Aside from the quantum statistical factors, this corresponds to the result given in |33j . (For 
a classical system obeying Maxwell-Boltzmann statistics one has to substitute n B {l + n B ) 
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by tim- For a fermionic system, e.g. made out of nucleons, obeying Fermi-Dirac statistics 
one would have instead n F {l — np) [0,00 

In our framework we are in the position to evaluate explicitly from the underlying mi- 
croscopic theory the correlation ((jF^(t)jF^(t'))) and then to compare with the intuitively 
suggested form of fl4.77| ). With the definition (|4.70|) of the fluctuating source one has for the 
correlator 



{{^{x,k,t)^(x',k',t'))) 



(4.78) 



dko dk' 
J 2tt 2tt 

k k o 



Ak k' A(k, k )A(k', k' )((AI(X, k, k Q )AI(X', k', k' ))) . 



(In our special case the spectral function A, and also I and f , do not depend explicitly on 
the center-of-mass variables as the background heat bath is assumed to be homogeneous 
and stationary.) The evaluation of the correlator ((A/A/)) is given in appendix |C] with the 
result ra 



((A/(X, k)AI(X', k'))) « (2tt) 4 5 (4) (X - X') [5 {4) {k + k') + 5^{k - k')] (l(k)f . (4.79) 

The sign is a reminder that one assumes the spectral function to be strongly peaked. 
Inserting ( |4.79[ ) in (|4.78|) one gets 



((^(f,fc,t)^(f,fc',t'))) = 

4(27r) 3 5 4 (X - X')5\k - k') (l(k,4)) 2 J ^ (A(k,k 



(4.80) 



for a quasi onshell spectral function A. Employing ( |4.39|) and assuming the width to be 
sufficiently small, the integral then can be evaluated as 



dko 
2^ 



(A(k, k 



Hence, one has 

((Ft(x , k, t)^(x', k', f ))) = S^(X - X') ((2tt) 3 ^ 3 ) (k - k')) I T 



(4.81) 



(4.82) 



S^(X-X') ((2n) 3 S {3 \k-k' 



1 (l(k,u§ 



In order to compare this result with the intuitive expression ( |4.77|) , one has to average over 
the space variables x and x! and identify the factor (27r) 3 <5 3 (&; — k') for k = k' with the 



volume of the system. In addition we employ the definition of the damping rate (|2.55| ), the 
property ( |2.59| ) and the generalized fluctuation-dissipation relation ( |2.63[ ). With some minor 
manipulations I? becomes 



It = 2T^(l + 2n B f 



2T [n B (l + n B ) + 



(4.83) 
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This corresponds nearly to the intuitively derived correlation strength of (|4.77|) , but deviates 



by an additive and constant number. The fluctuation (((A/e) 2 )) of the quantity around 
its average introduced in (|4.68|) thus takes the form 



(((A/ 5 ) 2 » = n B (l + n B ) + l -± n B (l + n B ) , (4.84) 

and therefore accounts not only for the thermal fluctuations (which would vanish for T — > 0). 
We are thus forced to conclude that f^(x, k, t) does not exactly correspond to the fluctuating 
onshell phase-space distribution! 

Before drawing any further conclusions, we also want to mention that one can also obtain 
the relation (|4.84 ) in a more direct way without going through the various approximations 



which lead from the exact equation of motion ([4.62 ) to the transport equation ( 4.66 ). One 



may define the fluctuating 'particle number density' of particles with momentum k as 



(n(*,t)> c :=± l<Pxf.(.r.k.l) (l.sr,) 



\ J >* J ^(2k ) (iD<(X,k) 



V 

o 

One can show that within the quasi-particle approximation, this definition is equivalent to 
the quantity (N(k,t))^ of ( [4.9|) . A lengthy calculation which is sketched in appendix |D| 
ends with the following expression ( |D12| ) for the average fluctuation of this quantity in the 



long-time limit: 

«(MM)) C ) 2 » = (n B (u$) 2 + n B (u,l)(l + n B (u>$) + - . (4.86) 



This result is equivalent to Q4.84Q 



We now have to point out what this lengthy enterprise in this section was good for. 
We have to admit that it has been our original intention to derive a Boltzmann-Langevin 
equation from first principles. In the beginning we have shown how to evaluate the particle 
density distribution and also how to evaluate its thermal fluctuation around its average. 
The latter has been obtained by means of a careful projection out of the equal-time four- 
point function. We were also tempted to introduce the quantity (N(k, t))^ in Q4.9Q which was 
expected to be the stochastic analogue to the average particle number. We then have derived 
(within the appropriate and necessary approximations) the kinetic transport equation for 
the phase-space distribution / describing the bosonic particles in the system. Its stochastic 
analogue then seemed to have the structure of the phenomenologically inspired Boltzmann- 
Langevin description. However, when calculating the average of the fluctuations of either the 
quantity /t or (N(k,t))^, we have now learned that their fluctuations are not only given by 
the expected thermal contribution, but instead are supplemented by an additional constant 
number. If the temperature is taken to zero (vacuum case) only this constant remains 
characterizing a kind of quantum fluctuation besides the thermal one. Hence, neither 
nor (N(k,t))^ do really have the physical meaning of what their labels have advertised. 
The appearance of additional quantum fluctuations can be readily traced back to the noise 
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correlator I which yields the Gaussian width of the noise fields £. Averaging over this 
Gaussian distribution also incorporates some quantum effects for the one-point functions 
as we have already seen when discussing the generalized (quantum) fluctuation-dissipation 
relation (|2.63| ). The coth term there (which is nothing but + 1) does not vanish for the 
vacuum case T = 0. Indeed, for low frequency modes, i.e. in the limit T ^> uj where ( |2.63| ) 
becomes identical to the classical fluctuation-dissipation relation ( p. 651 ) we ma Y neglect the 
constant 1/4 in ( l.84j ) compared to the large occupation number densities. For arbitrary 
frequencies, however, we have to keep this constant term. In addition we also want to note 
that the fluctuating source in the stochastic equation ( |4.69p is actually not distributed 
by a Gaussian process (as advertised above), as it is a quadratic quantity in the more 
'fundamental' (and Gaussian) noise £. 

We therefore conclude that there exist no real quantity which can account for a 'fluctu- 
ating particle number' on the level of a one-particle Green's function or two-point function. 
The information about the fluctuation of the particle number cannot be mimicked by any 
two-point function. It can only be obtained when considering and solving for the higher 
order Green's functions, e.g. the two-particle or four-point Green's function. On the other 
hand, this implies that a Boltzmann-Langevin equation as being advertised in the literature 
I2] , |54] ] can not be derived from any first principles, has not the meaning it was intro- 



duced for. This, however, does not necessarily mean that the phenomenological applications 
of Boltzmann-Langevin equations [|E| are not meaningful. It simply means that they stand 
for an intuitive advice incorporating some (thermal) fluctuations. 



V. REMARKS ON PINCH SINGULARITIES 

We have seen in the previous sections that modes or quasi-particles become thermally 
populated by a non-perturbative interplay between noise and dissipative terms entering the 
Kadanoff-Baym equations. It is non-perturbative as the equations of motion explicitly resum 
the self energy contributions by a Schwinger-Dyson equation defined on the real time path 
contour. In addition we have observed in eq. (|2.60|) that seemingly ill-defined expressions 
of the form 'O'/'O' are well-behaved in the sense of a weak coupling limit. Similar ill-defined 
expressions result from so called pinch singularities within the context of real time non- 
equilibrium field theory. It is the purpose of this section to demonstrate explicitly how 
pinch singularities are regulated within a non-perturbative context, though, at this point, 
we do not fully resolve the problem. Instead we try to give some (hopefully) helpful and 
intuitive remarks on this issue. The interpretation is rather similar to the observation seen 
in eq. ( |2.60p . We note that this section does not follow closely the mainstream of the other 



sections, but we believe it as being appropriate since our conclusions presented so far can 
shed also some light on this issue. 

A potential 'pinch singularity' arises (as an artifact) in strictly perturbative expressions 
of the form 

J 'dt 3 dUD I *(t 1 ,t 3 )A(t 3 ,U)D™(U,t 2 ) , (5.1) 
which for a (quasi) homogeneous and stationary system becomes by Fourier transformation 
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- D™\p,u)A(j>,uj)D™(p,uj) . (5.2) 

A describes some physical quantity (e.g. a self energy insertion) and is not further specified 
here. As D™ 1 contains a pole at u = E p — ie and Dq v contains a pole at u = E p + ie, 
the product of both in the above expression is ill-defined, if A(p, u = E p ) is not vanishing 
onshell. 

It was observed and proven by Landsman and van Weert that such ill-defined terms 
do cancel each other, to each order in perturbation theory, if the system stays at thermal 
equilibrium [[/J. Their arguments, however, rely solely on the KMS boundary conditions of 
the free propagators and self energy insertions, so that they do not apply for systems out of 
equilibrium. This severe problem arising for systems out of equilibrium was first raised by 
Altherr and Seibert |35[ . 



In a subsequent paper Altherr |36[ tried to 'cure' this problem by hand by introducing 
a finite width for the 'unperturbed' free CTP propagator D Q so that the expressions are at 
least defined in a mathematical sense. Such a procedure, of course, already represents some 
mixing of non-perturbative effects within the 'free' propagator. As we have seen from our 
discussions in the previous sections a damping term (resulting in a non vanishing width) 
and the associated noise guarantee that the propagator becomes thermally populated at 
long times. Hence, for any system which moves towards thermal equilibrium and thus 
behaves dissipatively, the full propagator must have some finite width (due to collisions or 
more generally due to damping). Plasmons behave as 'nonshell' modes |37|j38|] . Within 
his modified perturbative approach, Altherr |36[] also showed that seemingly higher order 
diagrams do contribute to a lower order in the coupling constant, as some of the higher 
order diagrams involving pinch terms will receive factors of the form 1/7", n > 1 reducing 
substantially the overall power in the coupling constant. (In his particular case Altherr 
investigated the dynamically generated effective mass within the standard 4 — theory. For 
the hard modes the onshell damping 7 is of the order o(g A T).) As a conclusion he raised that 
power counting arguments might in fact be much less trivial for systems out of equilibrium. 

Bedaque [|39| noted that pinch singularities are in fact an artifact of the boundaries 
chosen at to — > — 00. A dynamically evolving system is prepared at some finite time to. Time 
reversal symmetry is thus explicitly broken, so that the propagators have to depend on both 
time arguments explicitly before the system has reached a final equilibrium configuration. 
As long as the system is not in equilibrium (on a time scale of roughly I/7, 7 representing 
the relaxation rate), the propagator thus cannot be stationary. Furthermore he showed that 
by calculating the strictly perturbative expression (|5.1|) with free propagators starting from 



some finite to in the time integration, the occurring singularities in fact do cancel. As these 
explicit calculations are already rather involved, higher order diagrams in a perturbative 
series are expected to be even more complicated, and it is not clear whether such a procedure 
should really be necessary. 

In the following we want to emphasize on both approaches in more detail and try to 
elucidate on the problem with somewhat more physical intuition. In order to cure the 
pinch problem from a mathematical point, one has to resum the propagator to achieve a 
finite width. (Note again, that otherwise, in the long-time limit, the propagator would not 
become thermally populated.) Next we argue why Altherr has found that contributions of 
higher order have to contribute to lowest order. The basic reason is that any transport 
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process is non-perturbative and cannot be described by naive power counting. 

Let us start our considerations by looking at the Dyson-Schwinger equation of the form 



D = D + D ED 

= D + D ZD + A)£A)££o + • • • (5.3) 

on the real time path contour. By applying the so called Langreth-Wilkins rules, which are 



quite well known within the context of the Keldysh formalism ||40|| , one first finds for the 
retarded and advanced propagator 

pret = D rct + ^rct^rct^rct 

D rct 1 



av 




^av = £,av + ^av^av^av 



D av 1 
_^ D w = _ ^ / 5 4 n 

1-Z^ v £ av (Dg v )(- 1 )-S av v ; 

in symbolic notation. For D K we first stick to the perturbation expansion and do follow 
again strictly the Langreth-Wilkins rules 

D< = D< + D T et E ret D< + Dg Bt E < DJ Y + L^£ av .D av + 
D T Q et Z ret D T ct Z Ict D< + D5 8t E rot D5 Bt E < D5 Y + 

D<S av ^ v S av L) av + . . . . (5.5) 

Unavoidably, such a perturbative expansion would exhibit an infinite number of ill-defined 
expressions. (As a matter of fact, only the very first term in this expansion is regular, if 
E < (p, uj = E p ) is not vanishing onshell. Also terms invoking the product of DqD™* and 
DqDq v are as well ill-defined and lead to pinch singularities of type Q5.1|).) However, a 
simple rearrangement of the terms in the above series leads to 

D<=(l + L£ et £ ret + (DfS rct ) 2 + . . .) D< (l + S av D av + (S av D av ) 2 + . . .) 

+ Df (l + L>* et £ ret + (D* et £ ret ) 2 + . . .) £< (l + £> av £ av + (D av £ av ) 2 + . . .) £> av 

i i n ret n av 



1 - D5 ct S rct 1 - £> av £ av 1 - L>5 ct £ ret 1 - L> av £ av 

D ret ^ D rrt) (-!)£< p«r)(-l)^ ^av + ^retj^av _ ^ ^ 



The resummation of the series ( |5.5| ) of ill-defined terms results in a well-defined expression. 
This expression is very similar to Q2.49 ) which its meaning we already have discussed in 



detail in the previous sections. The first term stems from the possible boundary conditions 
specified at some finite time t = t and corresponds to the one in eq. (|2.49|) . As emphasized 
before, the second term guarantees that the modes become populated by the interaction (in 
our studied case by the interaction with the heat bath degrees of freedom) and thermally 
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populated in the long-time limit when the self energy parts S < ,S > do fulfill the KMS 
relation. (Note that the relations of this chapter are valid for arbitrary self energies. They 
need not be restricted to the case of a coupling to an external heat bath as we did in the 
previous sections.) 

Let us now concentrate more specifically on the last term in ( |5.6| ). If one assumes a 
homogeneous and stationary system and applies a Fourier transformation even for a non- 
equilibrium configuration (which is a contradiction, but let us stay to this assumption), one 
finds with T(p,u) : = ^(S> - £<) and ReE(p» = Re£ ret (p» = ReS av (p» 



p2 _ m 2 _ YleY, + iujY p 2 — m 2 — ReS — iuY 
1 



(p 2 — m 2 — ReE) 2 + uj 2 T 2 
= (_2t) ^ f — (5 7) 

In the limit of nearly vanishing width, i.e. assuming f, S > , S < to be proportional to some 
power g n in the coupling constant g, this results in 

J D ret S < J D av ^ -2msgn(uj)5(p 2 - m 2 ) lim ^ - > ^ < - < ^| . (5.8) 

As noted already in eq. (|2.60|) , this last expression stays well-behaved in the sense that the 
limit 

lim ( - S< - 1 (5.9) 

does exist and therefore is free of any pathological behavior, i.e. free of any pinch singularities. 
The physical meaning of the expression (compare ( [4.45| ) and discussion thereafter) 

becomes obvious (either in the limit g —> or not) when applying the equilibrium KMS 
conditions for the self energy parts 

nv(p,u) = I - > _ - < I — ► n B (uj). (5.11) 

This is just the Bose distribution function (see also eq. (|2.59p ). Hence, one has to interpret 
n T.iPi u ) as the 'occupation' factor demanded by the self energy parts. If the self energy 
parts do not fulfill the KMS relation, i.e. for a general non-equilibrium situation, this factor 
truly deviates from the Bose distribution. As derived in the previous section |TV| , the collision 
term ( [4.54J) in the transport equation (|4.58|) is proportional to the difference between the 
actual occupation probability n and the self energy occupation probability n^. 
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iE < (l + n) - iS>n = i(E< - E>)(n - n s ) = -2r(n - n s ) . (5.12) 

This term drives the system towards equilibrium and only there the two quantities n and 
ns become identical. 

For our arguments below, we want to note that the above expression (|5.6Q can also be 
cast in a different, but irritating form. We follow here the analysis of Baier et al (see 



also Carrington et al |42|), starting from the specification of a non-equilibrium distribution 
of the form 



Dq = —2irisgn(uj)5(p — m )n(p,uj) , (5.13) 

where h describes 'some' non-equilibrium distribution function. (We put 'some' in apostro- 
phes as it is not clear what kind of distribution one should really take!) Their result for 
doing the resummation ( |5.3| ) reads 



D< = -D iet n + D av h 

+ D TCt ((1 - n)S< + nS>) D av . (5.14) 



This expression is very similar to the one found (perturbatively) by Altherr and Seibert |3 
and Altherr ||36|| . It is easy to verify that (|5.14j) is identical to the second term in ( |5.6|) . (The 



first term in ( |5.6| ) obviously vanishes for an assumed homogeneous and stationary system as 
(-Dq *)*^ 1 ) = p 2 — m 2 and thus (p 2 — m 2 )5{p 2 — m 2 ) = 0.) For this to see we rewrite (|5.14| ) 
and employ 



D< = -nD TCt ((D™)^ - (D^y-V) D av 

+D rct S<D av + nD TCt ( £> - S <N ) D av 



-nD TCt (S rct - S av ) D av + nD rct (s> - E<) D av + D rct S<D av 
J D rct S < J D av , (5.15) 



as S ret - E av = £> - E<. 



The amusing thing to note is that in fact n drops out completely from ( |5.14j ) and therefore 



the form of the expression (|5.14j ) is in fact quite irritating. Hence, within this resummation, 



the propagator D < does not directly depend on the non-equilibrium distribution n. A 
possible dependence can only arise, if the self energy E < depends on h. So the question 
remains of how does h enter? 



Let us first try to understand the observations made by Altherr | 36fl . He has found, 
starting from some non-equilibrium distribution n, that higher order diagrams contribute to 
the same order in the coupling constant as the lowest order one. Indeed, in his investigation, 
the particular higher order diagrams where nothing but the perturbative contributions of the 
series in ( |5.5| ) for the dressed or resummed one-particle propagator D K . The only difference 
is that he has employed a 'free' propagator modified by some finite width in order that each of 
the terms in the series ( |5.5| ) becomes well defined. The reason for the higher order diagrams 
to contribute to the same order is now obvious. The initial out-of-equilibrium distribution n 
cannot stay constant during the evolution of the system as it has to evolve towards the Bose 
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distribution. If n — ns is of order o(l), it is obvious that there must exist contributions which 
perturbatively attribute to the temporal change of the distribution function and contribute 
to the same order as the one with employing ( |5.13| ). In fact, doing the resummation stated 



above as proposed by Altherr and also by Baier et al. ||41|| , we now see that the resummation 
would simply result in exchanging h by ns given by ( |5.10| ): 

n — > n^(p, uj) . (5.16) 

It should be clear by now that this 'change' stems in part from the temporal evolution of 
the distribution function. One can also interpret the substitution (|5.16|) as the change of 
the distribution dictated by the self energy operator (containing damping). 

However, this is still not the full solution: In a self interacting theory the self energy parts 
S < and S > do also evolve with time as they depend on the evolving distribution functions. 
Finally, in the long-time limit, they have to obey the KMS relation when the system has 
become completely equilibrated. Speaking more technically, the self energy operators have 
also to be evaluated consistently by the fully dressed and temporally evolving one-particle 
propagators. To make the argument more clear, we symbolically write 

n(t) « n B + (n(t ) - n B )e~ jt *- = ^> n B 

D(h,t 2 ) * 1 ^°° J D cq (t 1 -t 2 ) 

E(t u 1*)[D]*^V«(t l -t 3 )[D eq \ . (5.17) 

Within the time of about 1/F the functions under consideration do change in time and thus 
describe the non-equilibrium evolution (cf. ( |5.12| )). In particular, during the non-equilibrium 



stages, the one-particle propagator D as well as the self energy £ will depend on both time 
arguments explicitly. 

It is thus not clear what the assumption of quasi stationarity resulting in an evaluation 
of the contributions by Fourier technique really means. In fact, one could intuitively argue 
that if one can calculate all contributions with the simple Fourier description, one also would 
find higher order diagrams which do change the self energy operator. These would come 
from higher order corrections to the one-particle propagator employed inside the self energy 
operator. Then, the final answer for the 'fully' resummed propagator would be simply that 
one describing thermal equilibrium, i.e. 

^» - '- a ) (p2 _ m2 _^ q)2 + ^ ^M • (5.18) 

The basic reason for this is simply the fact that the system is out of equilibrium only 
for a finite interval in time. Applying standard Fourier techniques (by assuming quasi 
stationarity) simply washes out the intermediate non-equilibrium behavior. 

Our last arguments, though more intuitively, follow exactly the ones given already by 
Bedaque [|3S[| . The evolution of a non-equilibrium system is always in a strict sense non- 



perturbative. A full resolution of the pinch problem still has to be given. However, we 
feel that a full resolution should be handled more intuitively than the explicit (and still 
only strictly perturbative) calculations performed in []55| . The idea must be to write down 
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(within approximations) the full non-equilibrium equations for the one-particle propagators 
under consideration, i.e. the appropriate Kadanoff-Baym equations. The non-equilibrium self 
energy operators, which do enter there, must also be written down by means of the evolving 
dressed propagators. One carefully has to think about avoiding double counting, i.e. follow 
a skeleton expansion. Calculating some observable at some specified time depending on 
the propagators, one again has to stay in the appropriate order of the skeleton expansion 
to avoid double counting. Such a non-perturbative procedure should be free of any pinch 
problem right from the beginning. If one then succeeds in obtaining the correct kinetic 
and Markovian limit the actual non-equilibrium distribution(s) at that given time will enter 
explicitly p3 . 



VI. SUMMARY AND CONCLUSIONS 

We hope that our study has provided some new intuitive insight into (non-)equilibrium 
quantum field theory. In our discussions we have elucidated on the stochastic aspects inher- 
ent to the (non-) equilibrium quantum transport equations, the so called Kadanoff-Baym 
equations. For this we have started with a simple model, where we have coupled a free 
scalar boson quantum field to an external heat bath with some given temperature T. The 
coupling is achieved by some (linear) self energy operator which was assumed to be solely 
determined by the properties of the heat bath, i.e. it does not depend on the actual system 
variables. This operator can be thought as the effective interaction among the heat bath 
and the (now open) system when integrating out the heat bath degrees of freedom. We then 
showed that the emerging quantum transport equations for various two-point functions are 
just the familiar Kadanoff-Baym equations. The self energy operator was then divided into 
three distinct physical quantities ( |2.35| )-( p37| ) denoted by s, a and /. It turned out that the 
first term denoted by s is responsible for a dynamically generated mass shift, i.e. it renor- 
malizes the bare mass. The second term denoted by a is responsible for the damping of the 
propagating modes and is related via ( |2.55| ) to the usual plasma damping rate T. The last 
term denoted by I solely characterizes the (thermal and quantum) fluctuations inherent to 
the underlying transport process. For this to see we had rewritten the generating functional 
as a stochastic generating functional averaged over a (positively defined) weight of a Gaus- 
sian distributed noise field. This corresponds just to the influence functional approach of 
Feynman and Vernon. The emerging stochastic equations of motion can then be seen as 
generalized (quantum) Langevin processes. The generalized fluctuation-dissipation relation 
then states a direct relation between the damping rate T and the noise kernel /. Such a 
relation is of utmost microscopic importance as it dictates that the modes become thermally 
populated on the average. The process of thermalization results as an intimate interplay 
between dissipation and fluctuations. In turn, the Kadanoff-Baym equations are understood 
as an ensemble averaged description of the evolution of various two-point functions each of 
which belonging to a given noise sequence. 

What is changed, if we replace our toy model of a free system coupled to an external heat 
bath by a self-coupled and thus nonlinear closed system? In an interacting field theory of a 
closed system the Kadanoff-Baym equations formally have exactly the same structure as in 
our toy model. The important difference, however, is that the self energy operator is now 
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described fully (within the appropriate approximative scheme) by the system variables, i.e. it 
is expressed as a convolution of various two-point functions. Hence, an underlying stochastic 
process, as in our case an external stochastic Gaussian process, cannot really be extracted. 
However, we want to emphasize again that the emerging structure of the Kadanoff-Baym 
equations is identical. The carried out decomposition of the self energy operator into its 
three physical parts (mass shift s, damping a, and fluctuation kernel I) can immediately 
been taken over and thus keep their physical meaning also for a nonlinear closed system. At 
this point we want to mention two recent works in the literature: In a work by Csernai, Jeon 
and Kapusta \M\ it was shown that for a purely classical system one can also describe the 
microscopic equation of motions for each particle by means of an appropriately generalized 
Langevin equation, either in the case of an externally coupled heat bath or regarding a 
closed system. In another work, Crisanti and Marconi have taken the opposite direction of 
reasoning and have outlined how an effective action method valid also for non-equilibrium 
processes for describing the ensemble averaged properties of a system can be constructed if 
the underlying microscopic stochastic equation is of generalized (and nonlinear) Langevin 
type with white noise. The emerging equations of motion do resemble in their structure that 
of the real time description of (quantum) field theory. 

In section [TT1] we have emphasized that the long wavelength modes with momenta \k\ -C T 
and energies \u\ <C T do behave as classical propagating modes for a weakly interacting the- 
ory. Their evolution is described by means of a classical Langevin wave equation. This 
confirms the conclusion already stressed in [14] where the effective propagation of the soft 
modes within the self coupled </> 4 -theory had been investigated. The hard modes being (per- 
turbatively) integrated out do act as an external heat bath as they carry most of the energy 
of the system. The observation is that for a weakly interacting theory the average occu- 
pation number of the soft modes becomes large and approaches the classical equipartition 
limit. Their description can thus in fact be described classically supplemented by dissipative 
terms and fluctuating sources, which represent their (perturbative) interaction with the hard 
modes. Here we would like to mention that very recently there had been some interesting 
progress of how to describe the non-perturbative evolution of (super-)soft modes (on a scale 
of \p\ ~ g 2 T) in a non Abelian field theory by means of Langevin equations [|1^,0. We 
hope that our detailed analysis of the real time description of the soft modes within the 
Schwinger-Keldysh formalism can attribute to this subject. The understanding of the be- 
havior of the soft modes is crucial for the issue of anomalous baryon number violation in 
the hot electroweak theory. 

Semi-classical Langevin equations may not hold for a strongly interacting theory as for 
highly non trivial dispersion relations the frequencies of the long wavelength modes are 
not necessarily much smaller than the temperature. Still, when the soft modes become 
tremendously populated, as for example in the phenomena associated with the so called 
disoriented chiral condensates or during the inflationary epoch in an expanding universe, 
one again can argue that the long wavelength modes being coherently amplified do behave 
classically. In the work of Rajagopal and Wilczek ^7J, they had demonstrated that during 
the quench of a second order chiral phase transition the soft pion fields become highly 
populated and can be described by their classical equation of motions. On the other hand, 
these soft pion modes still do interact with the thermal environment which was not taken 
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into account (even classically) in [47 due to the finite resolution of the employed spatial 
lattice. To account for their interaction it has recently been proposed to supplement the 
classical equation by appropriate Langevin terms ||48|| . As the equations are clearly nonlinear 
and the initial configuration of the system is unstable, such a Langevin description can lead 



to various interesting phenomena for the possible evolution of the soft pion modes |H8 



In section [TV] we have derived the kinetic transport equation of the phase space dis- 
tribution function for the semi-classical particle regime. Such a transport description is 
strictly valid only within the quasi-particle picture (i.e. weak coupling regime) for a nearly 
homogeneous system. It has to be understood as an ensemble average over possible real- 
izations of the system. The self energy contribution / which characterizes the correlation 
kernel of the noise fields at the one-point level is already inherent in the Kadanoff-Baym 
equations which deal with the two-point functions. We have also studied extensively how 
to obtain the fluctuations of the phase space distribution. We tried to derive on the same 
footing a Boltzmann-Langevin equation from first principles from the underlying stochastic 
quantum transport equations. The obtained equation resembled in its structure the phe- 
nomenologically inspired Boltzmann-Langevin description which has a form of a standard 
transport equation supplemented by an additional randomly fluctuating term. A closer in- 
spection, however, showed that the inherent fluctuations do not solely correspond to the 
ones expected for a thermal system, but also contained an additional quantum correction. 
A microscopic justification of the phenomenological Boltzmann-Langevin approach can thus 
not be given. 

It has been the purpose of section [V| to demonstrate explicitly how so called pinch sin- 
gularities are regulated within the non-perturbative context of the thermalization process. 
Ill-defined pinch singularities do occur when higher order diagrams are evaluated perturba- 
tively by means of bare and unperturbed propagators, if the inserted occupation numbers 
are not the ones of thermal equilibrium. As we have outlined in section [V], the process of 
thermalization always comes along with dissipation so that the propagators naturally do 
acquire a finite damping width. This in turn regulates the pinch singularities. We have also 
commented in some detail on the existing literature on this topic. In particular we have 
argued that a strictly perturbative evaluation leads to some peculiarities in a naive count- 
ing scheme as raised by Altherr f36| . These are physically based on the non-perturbative 
character of thermalization. 



We have not touched in our study the important issue of decoherence (see e.g. [|49|,p0|l). 
Loosely spoken, decoherence of initial quantum mechanical correlations in a system occurs 
due to the stochastic forces of the (external or internal) heat bath acting randomly on 
the correlated trajectories and therefore destroying those correlations or phases with time. 
The density matrix describing the system becomes more and more diagonal with time as 
the off-diagonal entries are more and more washed out so that the individual trajectories 
will separate and behave nearly as classical. We want to refer the interested reader to 
the extensive literature on decoherence. The general principle can already be found in the 
seminal work of Caldeira and Leggett |T3, where it is shown how the density matrix of a 



quantum Brownian particle is decohered. This idea was recently further pursued by Elze |5T 



with special emphasize on the decoherence of partons in ultrarelativistic heavy ion collisions. 
There it was argued that the partons very fast lose their quantum mechanical 'affiliation' to 
the nucleons and then act as individual particles without any ongoing correlation to their 
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'mother' nucleoli. There also has been an extensive interest on the subject of decoherence 
in the evolution and decoupling of the 'wave function' of the universe [jl9l , p0|1 - The idea is 
that the whole universe acts as a huge (quantum) reservoir interacting with the individual 
subsystems. Accordingly, the individual evolution of the subsystems then decohere and lose 
their quantum mechanical correlations to the other subsystems. We also want to mention 
that the classical Einstein equation can effectively turn into 'Einstein-Langevin' equations 
if a quantum field is coupled to the metric and then integrated out. This might then lead 



to stochastic structure formation in the early universe [52 



We hope that the work presented here sheds some light on the relation between stochastic 
processes being described by equations of Langevin type and the underlying microscopic 
quantum theory. 
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APPENDIX A: REMINDER OF THE INFLUENCE FUNCTIONAL TECHNIQUE 

OF FEYNMAN AND VERNON 



A quantum mechanical system S (described by the variable x) interacts with a bath Q 
(described by the variable q). At some specified initial time t — t{ the combined system (S 
U Q) is described by the full density matrix 

Psuq(*<) = Psuqfa, ft, x'i, q[\ U). (Al) 

In the Schrodinger representation the (full) density matrix evolves in time according to 

PsuQ(xf,qf,x' f ,q' f ;t f ) 

= J dxidqidx'idq'i {l7(x f , q f ; x i} ft; t f , ^)psuq(^, ft; x[, g-; ti)f/ f (x-, ft'; x' f , q' f ; t f , U)} , (A2) 

where the evolution operator reads in the path integral representation 

x i if 

U(x f ,qf]Xi,qi]tf,U) = JVxJ Vq e* s ^ (A3) 

with 

*/ 

S[x, q} = J ds C(x(s), q{s)) = S s [x] + S Q [q] + S int [x, q] . (A4) 
u 

From the evolution equation ( |A2| ) for the (full) density matrix one is naturally lead to the 
"doubling" of the degrees of freedom in a real time description when calculating expectation 
values at some given later time t. In this sense the Keldysh path contour representation |3j 
follows straightforwardly (see below). 

The idea is now to formally integrate out the bath degrees of freedom q and thus to 
obtain an effective interaction S e s[x(s),x'(s)] describing the evolution of the system degrees 
of freedom x(s). 

Typically one furthermore assumes that the initial density matrix Psuq(^) is uncorrelated 
in its variables x and q, i.e. 

Psuq(^) = psfa, x' i: U) ® p Q (q i} q-; i f ), (A5) 

which can be motivated by assuming that the interaction £- mt (x, q) is adiabatically switched 
on at the time t{. (This assumption is not strictly necessary. However, it makes the following 
steps much more transparent. For a discussion of the treatment with inclusion of initial 
correlations between system and bath degrees of freedom see [pql .) In addition, the bath 
degrees of freedom are thought to stay in thermal equilibrium, i.e. Pq(^) = PQ ;e q(^), at 
some fixed temperature. 

Introducing the reduced density matrix p r as 

p r (x, x'; t) = Tr (g) {p(x, q; x\ q'; t)} = p s (x, x'\ t), (A6) 
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one finds that the expectation value of any operator A depending solely on the system 
degrees of freedom can be readily expressed in terms of p r : 

(A) = TT(x tq ) {A(x, x'\ t)p SuQ (x, q; x 1 , q'; t)} = Tr (a: ) {A(x, x'\ t)p r (x, x'\ t)} . (A7) 

Moreover one finds that its evolution in time can be put in the general form 

x f x 'f 

p r {x f , x' f ; t f ) = I dx,dx\ fVxf Vx' e i(5sW-5 s M) e is IF [*,*') pr ( Xh x '.. t .) ( A8 ) 



by introducing the influence functional Sw[x,x'] ]T7| 

e iSM*-*'] = Tl{q) {ug) q , tf} u)pQ{qu q >. tl )U { Q x,) \q' t , q f ; t, t,)} 

if if 

= J dqjdqidq'i J Vq J Vq' e ^S Q [ q ]-S Q [ q ']+S int [xM~S int [x', q ']) p ^ q ^ g ,. f ,y ^ 
n g'i 

Here x(s),x'(s) are treated as classical background fields. Thus one is led to say that p r 
evolves in time according to the effective 'influence' action 

S cS [x, x'] = S s [x] - S s [x'} + S w [x, x'}. (A10) 

It is clear that p r (i) evolves causally according to the history of the system and the bath. 

A notational simplification (the Keldysh representation 0) is achieved by introducing 
combined variables x c (s),q c (s) defined on the real time contour C (see the illustration in 
fig. 0) in the real time Green's function approach to finite temperature quantum field theory 
or non-equilibrium quantum field theory, i.e. 

/ x ( x(s) s G upper branch (1) / A11 \ 
xJs) = < } , \ i n i /„/ (All) 
L x (s) sG lower branch (2), 

and similarly for q(s), q(s'). The influence then takes the more compact form 

if 

elSnlxc] = J dqfdq . dq ' % J Vq J V q' e^/ CdSc { £ « (9c(s))+£int(:Cc(s) ' (?c(s)) }pQ(g i ,^;t i ), (A12) 

where the integration is defined on the time contour from t{ forward to t and back to ti. 
For the influence action Sip[x,x'] one finds the following general properties to hold [I7|: 



S JF [x,x']=-(S w [x',x})*, (A13) 
S lF [x,x}=0. (A14) 

Expanding S'if up to second order in x and x' (which corresponds to the regime of linear 
response), the general structure of Sif is given by p 
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t 

Sif[x,x'] ~ / ds F(s)(x(s) — x'(s)) 



i ; 

+ - J dsids 2 (x(si) - x'(si)) R(si, s 2 ) (x(s 2 ) + x (s 2 )) 
t 

+ I / ^ SlC ^ 2 ~ ( s i> ^2) (^(s 2 ) - z'(s 2 )) • (A15) 

Suppose now that £j nt takes the quite general form being linear in the system variables 

£int[x c , Qc) = x c (s)E (q c (s)) . (A16) 
Writing for Sjp[x) up to second order in x 

i „ r , | f C dsiF(si)a;c(si) — 5 f C dsids2X c (si)E(si,S2)aic(s2) , . . _ N 

exp-^i F [a; c ] w e R L J J (A17) 



it follows by comparing with ( |A12| ) 



'h\ 8e* SlF 
i J 8x c (si) 
'h\ 5 2 e% SlF 



F (si) = ( 7 ) 7T7-T = (S(ge(ai))>«, 



S(si, s 2 ) 



i / 5x c (si)5xi(s 2 ) 



= [<P(H(g e ( ai ))S(g c ( a2 )))> - (H(g e ( ai ))><S( 9c ( a2 ))>] (A18) 

where P means path-ordering along the contour C. F(s) can be interpreted as the average 
external force term due to the mean field generated by the average interaction of the bath 
variables Q with the system S. 

The path-ordering definition leads to the typical four self energy functions defined in real 
time: 

S 11 (ti,t 2 ) = -^((T(H(g(t 1 ))H(g(t 2 ))) - (mti)))(mt2))}) 

= e(ti - t 2 )s 21 (t 1; t 2 ) + e(t 2 - ti)s 12 (t l5 1 2 ) 

S 12 (ti,t 2 ) = -^((S(g(t 2 ))S(g(t 1 ))) - (Z(q(h)))(E( q (t 2 )))) = E<(t 1) t 2 ) 
Z 21 (t u t 2 ) = ~((Z( q (t l ))~( q (t 2 ))} - (E(q(t 1 ))} (E(q(t 2 ))}) = S>(t 1) t 2 ) 

S 22 (ti,t 2 ) = -^((T(S(g(t 1 ))H(g(t 2 ))) - (H(g(ti)))(H(g(t 2 )))) 

= 0{t x - t 2 )S 12 (t 1; t 2 ) + 6(t 2 - t x )E 21 (t 2 , t 2 ) (A19) 

Here T (T) stands for (anti)-time ordering. 
Defining the real Green's functions 
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Z R (tx,t 2 ) = E > (t 1) t 2 ) - E<(t 1 ,t 2 ) = -- ([E(q(tx)), E(q(t 2 ))]) 



S/(ti,t 2 ) = -(S > (t 1 ,t 2 ) + S<(t 1 ,t 2 )) 

= —(({mix)), E(q(t 2 ))}) - 2 (Sfofa))) (E(q(t 2 ))) ) , 



(A20) 



where E^(ti,t 2 ) = — S^(t 2 ,ti) and Ej(£i,t 2 ) = E/(/f: 2 ,ti), one finds after some algebra that 
in (| A 1 5|) i? and / are given by 



R(ti, t 2 



-^RihMWi-h) = -2Re(Z 11 (t 1} t 2 ))6(t 1 -t 2 ) 



I(h,t 2 ) = --Ztfuh) 



Im(S n (<i,i 2 )). 



(A21) 



From ( |A20| ) and ( |A21| ) one realizes that the function R is given by the (thermal) expectation 
value of the commutator of the E and corresponds to the usual (causal) linear response 
function. Physically it incorporates the dissipative effects between the bath and system 
degrees of freedom. 

On the other hand the function / which specifies the imaginary part in the influence 
action ( |A15|) is given by the (thermal) expectation value of the anticommutator (i.e. the 
symmetrized correlation function) and thus corresponds to a measure of the fluctuations in 



'5 — (S)) of the interactions stemming from the bath variables |L1]]. These act as random 



"kicks" on the actual trajectory and can be interpreted as a resulting stochastic force [JT8 
To see this in more detail, one defines the real stochastic influence action 



Sif[x, x', £] = Re(SiF[x, x']) + / ds £(s) (x(s) — x'(s)) 



(A22) 



where £(s) is interpreted as an external force, randomly distributed by a Gaussian distribu- 
tion with zero average: 



P[£00] = ^exp (-^j ds 1 ds 2 asx)r 1 (s 1 ,s 2 )as2) ) j 



(A23) 



The influence functional Sw [x, x'} is regained as the characteristic functional over the average 
of the random forces. 

e iSM*,x'] = _ = J Di P[i{s)} e^iFK*',?]. (A24) 

The imaginary part I then defines the correlation of the random forces 

(asins2))^ = ni( Sl ,s 2 ). (A25) 
Combining ( |A10|) , (|A15|) and ( |A22|) , the real effective and stochastic action 
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t 

S c s[x,x',^} ~ Ss[x] — Ss[x'} + J ds F(s)(x(s) —x'(s)) 

u 

i ; 

+ - J dsids 2 (x(si) - x'(si)) R(sx, s 2 ) (x(s 2 ) + x'(s 2 )) 
u 

t 

+ J ds £(s) {x(s) - x'{s)) (A26) 

ti 

governs the effective propagation of the system. 

In the case of a single particle this effective action will result in a generally retarded 
(non Markovian), stochastic Schrodinger equation for the propagation of the (stochastic) 
wave function []54[ . In the quasi-classical regime (or approximation) an effective Langevin 
like equation 



SS x [xci\ 



- F(t) - J ds'R(t, t')x d (t') = ^(t) (A27) 



5x cl (t) 

describes accordingly the evolution x c i(t) of the classical particle. (Classical) Brownian 



Motion |18| is recovered if the spectral function of the system allows for approximating the 



kernels R and / as 



R(s, s') = —m-fd^s, s') , 

M(s, s') = 2m 1 kT5(s, s') . (A28) 



APPENDIX B: GRADIENT EXPANSION 

Here we prove the following relation for the Wigner transform of a convolution integral: 
JdSi e iku Jd% S(X + u/2, y) G{y, X - u/2) = e- i0 S(X, k) G{X, k) (Bl) 
with G and S being the Wigner transforms of G and E, respectively. The operator O is 



defined as JJH 



2 

We start with the Fourier transform of the r.h.s. of (|B~T|): 

d\ 



<>--=\ {d^ - dM) ■ (B2) 



J^e-^e^(X,k)G(X,k) 

e -iku e -io U e ikv E ( X + v / 2j x- v/2) 



(2rcY 
d 4 k 



{2<kY 
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x fd\v e ikw G(X + w/2,X - w/2) 

J -0- e -iku J fa fa e ikv-lvd° e ihw+^ E ^ x + V / 2 ,X- v / 2 ) 

xG(X + w/2,X - w/2) 

k ~ iku f d% d\v e ikv e ikw Y,{X + (v + w) /2, X - (v - w) /2) 



(27T) 4 

xG(X + {w- v)/2, X-{w + v)/2) 
= Jd% E(X + u/2, X-v + w/2) G(X -v + w/2, X - w/2) 

= JdSj E(X + w/2, y) G(y, X - u/2) . (B3) 
Inverting the Fourier transform with respect to k yields (|HT|). 



APPENDIX C: THE CORRELATOR {{AI(X, K)AI(X', K 1 ))) 

In this appendix we evaluate the (ensemble averaged) correlator ((Al(X, k)Al(X', k'))), 
which enters the calculation of the correlator ((^(x, k, i)^(x', k', t'))) (eq. ([4.78Q ) of the 
fluctuating source !F^(x, k,t) in the kinetic stochastic equation ( |4.69| ). 

From the definition of the Wigner transformation Q4.25] ) we have 

((AI(X,k)AI(X',k'))) = 

J d'y d 4 y> e *»e*V ((AI(X + y/2, X - y/2)AI(X' + y'/2, X' - y'/2))) . (CI) 

Inserting the definition ( |4.63j ) of AI and making use of the Gaussian property of the noise 
distribution, i.e. the decomposition of the four-point correlator 

(ii(xi)i(x 2 )i(x 3 )i(x,))) 
= ((£(xi)£(x2)))((£(x 3 )£(x 4 ))) + ((axMxsMiixMx,))) + ((iixxnx^uiix^ixs))) 

= I(x 1 ,x 2 )I(x 3 ,x 4 ) + I(x 1 ,x 3 )I(x 2 ,x^) + I(x 1 ,x 4 )I(x 2 ,x 3 ) (C2) 
one finds 

((AI(X + y/2, X - y/2)AI{X' + y'/2, X' - y'/2))) 
= I(X + y/2, X' + y'/2)I(X - y/2, X' - y'/2) 

+ I(X + y/2, X' - y'/2)I(X - y/2, X' + y'/2) 
= I(X -X' + (y- y')/2)I(X - X' - (y - y')/2) 

+ /(X -X' + (y + y')/2)I(X - X' - (y + y')/2) . (C3) 

In the last equation we have used the fact that in our special case discussed in the main 
chapters the noise kernel / like all our self energies only depends on the difference of its space- 
time variables as the background heat bath is assumed to be homogeneous and stationary. 
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One now has to express the noise kernels / by their Fourier transforms and evaluate the 
spatial integrals over y and y'. One then gets 



((AI(X,k)AI(X',k'))) = 
d 4 p I(k + p/2) I{k - p/2) \5 (4) (k + k') + 5 (4) (k - k') 



(C4) 



-ip(X-X') 



As a further step, we now want to consider the Markovian limit of this expression, so that 
within this approximation the correlator will become effectively white and local. For this 
we consider the expansion 



I(k + p/2)=I(k) + 

I(k + p/2)=I(k)-^.p(k) 



+ o{p 2 ) 



(C5) 



fc=0 



+ o{p 2 ) 



k=0 



and truncate it after the linear term. Such an expansion should be acceptable if I(k) is a 
smooth function in k. Inserting ( pop into ( p4]) , one notices that the linear corrections do 
in fact drop out. Neglecting all quadratic and higher corrections, one thus receives for the 
correlator the following approximate form 

((AI(X, k)Al(X', k'))) « (2tt) 4 5 (4) (X - X') [8^(k + k') + 5 {A \k - k')] (l(k)) 2 . (C6) 

As a final remark we want to mention that the expansion (05) and hence the approximation 
for obtaining ( |C6|) is consistent with the spirit of the gradient expansion carried out in 
section [IV] for obtaining the quasi-classical kinetic transport equation: If the (phase space) 
distributions are sufficiently smooth, i.e. AXAk ^> 1, the terms neglected in are of 
order ~ o(l/(AXAA;) 2 ) and thus small. 



APPENDIX D: THE ENSEMBLE AVERAGED FLUCTUATION OF (N(K,T))^ 

In section [TV] we have defined the fluctuating 'particle number density' of particles with 
momentum k as 

1 



(n(k,t))r.= y I d s xfs(x,k,t) (Dl) 



vl d3x I^ {2ko) ( i5 f( x >v 



v 7^ (2A;o) / du ° elkouo I £x I d3 y e ~ rk " e 



x (iDf (t + M /2, x; t - u /2, y}) . 



In the long-time limit and within the quasi-particle approximation it is easy to see that the 
ensemble average of this quantity just gives the onshell Bose distribution function, i.e. with 
( p.56| ), ( p.59| ) and (|2.61|) one immediately recovers 
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«MM)>€» 



dk _ 



(2fc ) (iD < (k, ho) 



(D2) 



For its fluctuation we have to evaluate the following expression: 

(«n(£,f)) e <n(£,f)>«» 



(D3) 



l )TH I W^( 2k o) I du e 



2tt 



dun e 



ik u 



d 3 x / d^e-^e^ / d 3 x / dPye^e** 5 



x «£><(* + M /2, f ; * - mo/2, f) (* + «oA £; * - «o/2, y))) . 

Employing now the solution ( [4.11 ) in the long-time limit for the fluctuating two-point func- 
tion Df, and also using ( |C2|) , one has 



«D<(1,2)D<(3,4)» = D<{1,2)D<{3,4) 

r D ret IDa ,y\ 3) 

D TCt ID av ) (1,4) 



(D4) 



D TCt ID™ 
D rct ID av 



(2,4) 
(2,3) 



One can express the last two terms by means of the two-point functions D K and D > defined 
in ( p. 19 ), (j2.2Q|) . For this we have to write down the Kadanoff-Baym equation for D > being 
the analogue to ([2.48 ), which has not been mentioned in the main text so far: 



(-□ - m 2 - s - o)D > + (-a + U)D a 
Hence, in the long-time limit, one has 



0. 



'-i) D TCt ID &v 



1 



(zr 



D- 



and thus 



(( J D<(1,2) J D<(3,4)))= D<(1,2)D<(3,4) 



(D5) 



(D6) 



(D7) 



+ i [(£>< + D>)(1, 3)] [(D< + D>) (2, 4)] 

+ i[(D< + D>)(l,4)][p< + J D>)(2,3)] . 

When inserting ([D7D into Q , one recognizes that the first term on the rhs. of (|D7| ) just 
yields 

«MM)) f » 2 » (r^)) 2 . 

For the other two terms we know that in the long-time limit the (average) propagators 
become homogeneous and stationary, i.e. D(x, t x ; y, t y ) — D(\x — y\,t x — t y ). Evaluating all 
spatial integrals, one is left with 
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«<n(M)>e<n(M)>€» - (((n(k,t))sf 



(D8) 



:-i) 



dkn 
27 



(D < + D > ) (k, 



du e ikouo 
u + u 



dkn ~ 



2tt 



k / du e ifco "° 



(D < + D > ) (k, 



Uo + Uq, 



When expressing the propagators in time by their Fourier transforms, one finds after some 
further manipulations 



<«n(M)M»(M)>*» - «WM))*» S 



dkn 



2tt 



du 



2)/^(M 2 7T- (D < + D > )(k,u) (D< + D > )(k 1 uj-2k 



2vr L 



In the quasi-particle onshell limit, we recall 

D < (k,iu) = D > (-k, -uj) = D > (k,-uj) 



and thus 

D < (k 1 u) + D > {k ) u) 



-— (1 + 2n fl (u$) (5{u - W 2) + 5(a; + w 



Inserting this in ( |D9| ) , one receives the final expression 



«<n(M)>*<n(M)>*» - (((n(k,t))^y » - 1 + 2n B (u; 



n B (a# (l + + J . 



(D9) 



(D10) 



(Dll) 



(D12) 
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FIGURES 
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FIG. 1. Schematic picture for the interaction of the system with the heat bath. 
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FIG. 2. Schwinger-Keldysh time contour path C for the variables x(s) and x'(s) running from 
s = U to +00 and back to h. 
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